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Progress  report  from  first  year  of  JSEP  fellowship  (October  1997  to  present) 

The  last  year,  aided  by  the  Graduate  Fellowship  I  received  from  the  Joint 
Services  Electronics  Program,  has  been  a  fruitful  one.  In  the  following,  I 
summarize  the  progress  I  have  made  in  my  education  and  research. 

In  the  past  year,  my  research  has  focused  on  Bose-Einstein  condensates  stored 
in  an  optical  trap.  Briefly,  Bose-Einstein  condensation  is  a  low-temperature 
quantum  statistical  phenomenon  that  occurs  in  gases  of  bosons,  which  are 
particles  which  are  allowed  to  coexist  in  a  quantum  state  (unlike  fermions, 
such  as  electrons).  Once  the  temperature  of  the  gas  is  low  enough,  the 
wavefunctions  of  the  particles  comprising  the  gas  begin  to  overlap,  and 
further  cooling  places  nearly  all  the  particles  in  the  gas  into  a  single,  giant 
matter  wave.  Last  year,  our  group  made  a  significant  advance  in  the  field, 
trapping  a  Bose-Einstein  condensate  in  the  focus  of  an  infrared  laser  beam. 

This  trap  has  allowed  us  to  conduct  four  major  experiments  in  the  past  year. 
First,  we  were  able  to  place  an  optically  trapped  Bose-Einstein  condensate  in  a 
large  magnetic  field  and  observe  a  Feshbach  resonance  which  drastically 
modifies  the  interactions  between  atoms  in  our  Bose  condensates  of  atomic 
sodium.  We  were  able  to  show  that  this  resonance  significantly  modifies  the 
properties  of  the  Bose-Einstein  condensates  which  we  are  studying,  a  tool 
which  in  the  future  should  allow  scientists  to  arbitarily  "tune"  their  Bose 
condensed  sample  [S.  Inouye,  et  al..  Nature  392,  151  (1998)]. 

Second,  creating  traps  using  laser  light  has  the  benefit  of  being  able  to  define 
very  narrow  features  in  the  container  which  holds  the  Bose  condensate.  We 
used  this  ability  to  demonstrate  that  Bose-Einstein  condensates  could  be 
produced  in  a  novel  way  —  just  by  compressing  a  trapped  gas  above  the  Bose- 
Einstein  condensation  transition  temperature  using  a  few  milliwatts  of 
infrared  light.  This  project  was  especially  satisfying  as  an  elegant  case  study  of 
thermodynamics,  illuminating  basic  concepts  such  as  reversibility,  entropy, 
and  phase-space  density  [D.M.  Stamper-Kurn  et  al.,  Phys.  Rev.  Lett,  (in  press)]. 

Third,  we  used  the  optical  trap  to  confine  Bose-Einstein  condensates  that  were 
transfered,  using  magnetic  resonance  techniques,  into  combinations  of 
internal  hyperfine  states.  Placing  the  trapped  atoms  into  several  atomic  states 
leads  to  a  new  kind  of  quantum  fluid,  a  three-component  spinor  Bose 
condensate.  As  a  first  exploration  of  this  new  quantum  fluid,  we  studied  the 
formation  of  spin-domains  in  the  ground  state  of  spinor  condensates,  varying 
the  magnetic  fields  in  which  they  reside  [J.  Stenger  et  al..  Nature  (submitted)]. 

Fourth,  we  are  currently  exploring  Bose  condensates  composed  of  just  two 
components,  where  we  have  discovered  the  spontaneous  formation  of 
metastable  states  in  which  the  Bose  condensate  settles  into  a  high-energy 
state,  and  cannot  relax  to  the  true  ground  state  of  the  system  for  many 
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seconds.  We  are  investigating  the  mechnisms  by  which  this  metastable  state 
decays  to  the  ground  state. 


Finally,  in  the  past  year  I  collaborated  with  Martin  Naraschewski,  a  researcher 
at  the  Institute  for  Theoretical  Atomic  and  Molecular  Physics  (ITAMP)  at 
Harvard  University,  to  write  a  theoretical  paper  about  the  finite-temperature 
behaviour  of  trapped  Bose  gases.  The  paper  gives  analytical  expressions  for 
the  density  profiles  and  condensate  fractions  of  such  gases  which  should  be  of 
immediate  use  to  experimentalists  in  the  field  [M.  Naraschewski  and  D.M. 
Stamper-Kurn,  Phys.  Rev.  A  (in  press)]. 

In  addition,  I  have  had  the  opportunity  to  present  my  work  and  that  of  my 
colleagues.  In  July,  I  traveled  to  Paris,  France  where  I  visited  the  laboratories 
at  the  Ecole  Normale  Superieure,  a  strong  center  of  atomic  physics.  There,  I 
presented  recent  work  during  an  hour-long  seminar.  Recently,  I  presented  a 
poster  at  the  International  Conference  of  Atomic  Physics  (ICAP)  in  Windsor, 
Canada  [Poster  B24,  ICAP16  Book  of  Abstracts,  Windsor,  Canada  (1998)]. 

Last,  I  had  two  opportunities  to  advance  my  education  as  a  scientist  and 
future  educator.  This  spring  I  hired  an  undergraduate  student,  Jeff  C.  Gore,  to 
join  our  laboratory  for  the  summer  and  the  coming  year.  This  gave  me  the 
opportunity  to  meet  with  many  students,  and  to  learn  how  to  mentor 
newcomers  in  the  field,  an  experience  that  will  be  valuable  in  my  future 
career.  Also,  this  summer,  I  attended  the  "Enrico  Fermi"  Summer  School  on 
Bose-Einstein  condensation  in  Varenna,  Italy,  where  I  learned  more  about  the 
field,  and  had  the  chance  to  meet  and  work  with  my  colleagues  from  around 
the  world. 

In  conclusion,  it  has  been  a  very  busy  and  productive  year.  I  thank  the  JSEP 
Graduate  Fellowship  Program  for  making  this  work  possible. 

Yours, 

Dan  Stamper-Kurn 
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OPTICALLY  CONFINED  BOSE-EINSTEIN  CONDENSATES 


J.  Stenger,  D.M.  Stamper-Kurn,  M,R.  Andrews,  A.P.Chikkatur,  S.  Inouye,  H.-J. 

Miesner,  and  W.  Ketterle 

Daparfment  of  Physics  and  Research  Laboratory  of  Electronics, 

Massachusetts  Institute  of  Technology,  Cambridge,  MA  02139,  USA 

to  appear  in  JLTP 
12.  August  1998 


With  an  optical  dipole  trap  it  is  possible  to  confine  Dosc-Einstein  condensates  in  different  hyperfinc 
states  and  in  arbitrary  magnetic  bias  fields,  thus  overcoming  two  major  limitations  of  magnetic 
traps.  In  this  rewier  paper  we  characterize  the  properties  of  such  a  dipole  trap  and  we  summarize 
experiments  which  made  use  of  the  new  experime  ntal  i)ossibilities,  including  the  re  versible  formation 
o/  a  Detse -Einstein  condensate,  the  observation  of  Feshbach  resonances  in  sodium  and  the  ground 
state  properties  of  spinor  Be>sc-Einstein  condensates.  Finally,  we  present  some  new  rx suits  on  the 
shape  of  magnetically  trapped  and  ballisiically  expanding  condensates. 

PACS  numbers:  03.75  Fi.  05.30  .Jp,  32.80  Pj,  64.60 


1.  INTRODUCTION 

Three  years  after  the  first  observations  of  Bose-Eiustein  condensation  (DEC)  in  dilute  atomic 
gases^*'^  activities  in  this  new  field  are  still  increasing,  partially  due  to  the  many  new  experiments 
which  have  been  able  to  generate  Bose-Einstein  condensates.  All  studies  of  this  phenomenon  have 
been  performed  in  magnetic  traps.  These  traps,  in  combination  with  rf-induced  evaporation,  appear 
ideal  for  cooling  and  trapping  atoms  at  very  low  temperature  and  might  develop  into  workhorses 
for  nanokelvin  atomic  physics. 

However,  magnetic  traps  have  severe  limitations.  They  require  large  scale  inhomogeneous 
magnetic  fields  which  might  interfere  with  applications  in  precision  atom  optics.  For  example, 
in  the  first  demonstration  of  an  atom  la.ser,*  coherent  atomic  pulses  were  coupled  out  into  an 
inhomogeneous  magnetic  field,  which  served  to  confine  the  remaining  condensate.  Thus,  during 
propagation,  the  pulses  were  exposed  to  Zeeman  shifts.  While  these  shifts  were  mitigated  by 
producing  mp  =  0  atoms,  quadratic  Zeeman  shifts  may  preclude  precision  experiments  on  such 
pulses. 

The  unprecedented  control  over  the  motion  and  position  of  neutral  atoms  realized  by  BEC 
should  offer  major  advantages  for  metrology.  However,  magnetic  traps  cannot  trap  atoms  in  the 

=  0  .state,  which  are  preferable  for  atomic  clocks  and  other  precision  experiments,  thus  limiting 
the  use  of  trapped  condensates  for  metrology. 

More  generally,  magnetic  traps  can  confine  only  weak-field  seeking  hyperfine  states.  Since 
the  atomic  ground  state  is  always  strong-field  seeking,  weak-field  seeking  states  can  inelastically 
scatter  into  the  ground  state  (dipolar  relaxation)  resulting  in  heating  and  trap  loss.  This  restriction 
to  weak-field  seeking  states  also  limits  the  use  of  BEC  for  studies  of  various  atomic  properties 
which  depend  crucially  on  the  hyperfine  state,  such  as  collision  resonances,  or  for  examinations  of 
condensates  composed  of  various  combinations  of  hyperfine  states. 


Optically  confined  Bose  Einstein  condensates 

These  limitations  have  led  us  to  the  development  of  an  all-optical  trap  for  Bose-Einstein 
condensates.  In  this  review  paper,  we  summarize  the  basic  properties  of  this  trap  (Section  2).  and 
then  illustrate  various  new  possibilities  for  experiments  on  Bose-Einstein  condensation,  such  as  the 
following; 

•  Optical  traps  allow  precise  spatial  (micrometer)  and  temporal  (microseconds)  manipulation 
of  Bose-Einstein  condensates.  This  should  allow  the  realization  of  box  traps,  atom  guides, 
and  optical  lattices.  We  have  used  the  spatial  resolution  afforded  by  optical  traps  to  create  a 
new  “dimple'’'  trap  in  which  BEC  was  created  adiabatically  and  thus  reversibly  (Section  3). 

•  Optical  traps  have  a  new  external  degree  of  freedom:  they  can  be  operated  at  arbitrary 
external  magnetic  fields.  We  have  used  this  feature  for  the  observation  of  Feshbach  resonances 
for  strong  field  seeking  states  of  sodium  which  cannot  be  confined  magnetically  (Section  4). 

•  Optical  traps  offer  a  new  internal  degree  of  freedom:  the  orientation  of  the  atom's  magnetic 
moment.  This  resulted  in  the  generation  of  spinor  condensates,  condensates  which  populate  all 
three  hyperfine  states  of  the  F  =  1  state  of  sodium  (Section  o)  and  possess  a  three-component 
vectorial  order  parameter. 

Finally,  in  Section  6.  we  present  some  unpublished  quantitative  results  on  the  shape  of  mag¬ 
netically  trapped  Bose-Einstein  condensates,  both  in  the  trap  and  in  ballistic  expansion. 

2.  OPTICAL  CONFINEMENT  OF  A  BOSE-EINSTEIN  CONDENSATE 

Bose-Einstein  condensation  in  an  optical  trap  has  been  the  ultimate  goal  of  optical  cooling 
and  trapping  techniques  and  has  provided  one  motivation  for  the  development  of  sub-recoil  cool¬ 
ing  techniques, the  development  of  various  optical  dipole  traps, and  for  pursuing  Raman 
cooling^^'^'^  and  evaporative  cooling^  *  in  such  traps.  The  highest  phase  space  density  achieved  by 
purely  optical  means  was  a  factor  of  400  below  that  required  for  Bose-Einstein  condensation.^^ 

Rather  than  exclusively  using  optical  means  to  create  BEC,  we  have  adopted  a  different  ap¬ 
proach  for  studying  optically  trapped  condensates:^^  first  atoms  were  evaporatively  cooled  in  a 
magnetic  trap,  and  then  transfered  into  an  optical  trap.  This  approach  circumvents  many  difficul¬ 
ties  usually  encountered  with  optical  dipole  traps,  such  as  trap  loss  from  heating  processes  which 
are  proportional  to  laser  power  (off-resonant  Rayleigh  scattering),  or  which  lead  to  an  exponential 
heating  with  an  inverse  time  constant  proportional  to  the  laser  power  (heating  due  to  fiuctuations 
in  the  intensity  and  position  of  the  laser  beam^^).  Since  the  temperature  of  atoms  is  reduced 
through  rf-evaporation  by  a  factor  of  100,  only  milliwatts  of  laser  power  are  needed  as  compared 
with  several  watts  used  to  directly  trap  la.ser-cooled  atoms.  Furthermore,  since  the  cloud  shrinks 
while  being  cooled  in  the  magnetic  trap,  the  transfer  efficiency  into  the  small  trapping  volume  of 
an  optical  dipole  trap  is  increased. 

The  first  step  in  this  approach,  the  creation  of  a  Bose-Einstein  condensate  in  a  magnetic 
trap,  was  similar  to  our  previous  work.^^  Sodium  atoms  were  optically  cooled  and  trapped,  and 
transferred  into  a  magnetic  trap  where  they  were  further  cooled  by  rf-induced  evaporation.  This 
resulted  in  condensates  of  about  20  x  10^  atoms  in  an  elongated  cigar-shape  trap  due  to  our 
magnetic  trapping  configuration. 

The  optical  trap  was  formed  by  focusing  a  single  near-infrared  la.ser  beam  into  the  center  of 
the  magnetic  trap  along  the  axial  direction.  This  realized  the  simple  single-beam  arrangement  for 
an  optical  dipole  trap.'”^® 

The  parameters  of  an  optical  dipole  trap  formed  at  the  Gaussian  focus  of  a  single  red-detuned 
la.ser  beam  are  characterized  by  the  wavelength  A  of  the  light  (and  thereby  the  detuning  from 
resonance  5),  the  total  la.ser  power  F,  and  the  beam-waist  radius  u^o-  The  trap  depth  is  proportional 
to  P/w’^S.  Since  trap  depths  on  the  order  of  a  microkelvin  are  sufficient  to  trap  sodium  condensates, 
we  could  use  infrared  light  (A  =  985  nm)  which  was  far-detuned  from  the  sodium  resonance  lines 
(at  589  nm).  At  this  wavelength,  using  a  typical  beam  radius  of  wq  =  6//m,  the  trap  depth  is 
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1  f.iK/m\\\^^  and  the  aspect  ratio  of  the  atom  cloud  is  27.  Due  to  the  large  detuning  of  the 
infrared  beam,  the  spontaneous  scattering  rate  is  small  (5  x  per  //K  trap  depth).  None 

of  our  experiments  so  far  has  been  limited  by  the  spontaneous  scattering  rate  -  therefore  there 
is  currently  no  strong  incentive  to  reduce  this  rate  even  further  by  using  far-infrared  traps^^  or 
blue-detuned  traps.^^ 

Condensates  were  transferred  into  the  optical  trap  by  holding  them  in  a  steady  magnetic  trap 
while  ramping  up  the  infrared  laser  power,  and  then  suddenly  switching  off  the  magnetic  trap. 
Nearly  complete  transfers  of  the  atoms  were  observed,  with  condensate  numbers  in  the  optical  trap 
as  high  as  10'  atoms.  For  various  experiments,  we  have  produced  pure  optical  traps  with  beam 
radii  ranging  between  about  6  and  20  fim.  and  with  infrared  powers  ranging  from  about  5  to  50 
m\V. 

Condensates  in  the  optical  trap  are  usually  pure,  i.e.  the  condensate  fraction  is  very  high. 
This  is  because  the  optical  trap  has  such  a  small  volume  that  the  number  of  uncondensed  atoms  is 
very  low.  An  estimate  for  this  number  can  be  obtained  by  assuming  that  the  temperature  T  of  the 
cloud  is  about  1/lOth  of  the  trap  depth  (the  optical  trap  has  ‘'built-iiv*  evaporative  cooling  due  to 
the  limited  trap  depth).  The  number  of  non-condensed  atoms  is  approximately  [ksT/h^')^,  which 
for  typical  conditions  is  a  few  times  10*  atoms,  in  quantitative  agreement  with  our  measurements. 
With  a  condensate  of  about  5  million  atoms,  this  estimate  indicates  a  condensate  fraction  of  more 
than  99  %.  “Pure”  condensates  in  magnetic  traps  usually  have  a  larger  non-condensed  fraction 
due  to  the  smaller  trapping  frequencies  and  also  due  to  the  technical  difficulty  of  adjusting  a 
small  trap  depth  using  rf-induced  evaporation  since  the  trap  depth  is  sensitive  to  stray  magnetic 
fields.  In  the  future,  the  condensate  “purification”  easily  obtained  with  optical  traps  may  be  used 
in  experiments  to  observe  various  predicted  zero-temperature  phenomena,  such  as  collapses  and 
revivals  of  collective  excitations,  small  shifts  of  excitation  frequencies  due  to  quantum  depletion 
which  might  otherwise  be  masked  by  finite-temperature  effects,  or  generally  for  high-coherence 
atom  lasers. 

The  lifetime  of  optically  trapped  condensates  was  studied  by  measuring  the  number  of  con¬ 
densed  atoms  in  time-of-flight  images  after  a  variable  storage  time  in  the  optical  trap.  The  observed 
loss  rates  per  atom  in  the  optical  trap  ranged  from  4  s“*  at  a  peak  density  Uq  =  3  x  10*^ cm”^  to 
less  than  1/10  s”*  at  7?o  =  3  x  10**cm“'^.  The  obser\'ed  decay  was  predominantly  due  to  three- 
body  recombination.  The  three-body  loss  rate  constant  was  determined  for  sodium  condensates  as 
A3  =  1.1(3)  x  10"^^cm^s"*.  Fluctuations  of  the  optical  trapping  potential  did  not  severely  limit 
the  lifetime  of  the  samples,  which  at  lower  densities  exceeded  10  s. 

From  these  lifetime  mesurements,  we  can  also  e.stimate  an  upper  limit  for  the  heating  caused  by 
beam  jitter  and  spontaneous  scattering  in  the  optical  trap.  By  fitting  the  observed  lifetime  curves, 
we  estimated  a  density  independent  loss  rate  of  0.03(2)  s""*.  Assuming  a  trap  depth  of  about  5 
//K.  this  implies  a  heating  rate  of  about  150  nK/s.  which  is  larger  than  the  25  iiK/s  heating  rate 
one  would  expect  from  spontaneous  scattering  alone.  Another  estimate  of  this  heating  rate  was 
obtained  by  shining  the  infrared  light  onto  a  sample  of  gas  in  the  magnetic  trap,  and  monitoring 
the  temperature  of  the  gas  in  this  combined  potential.  Accounting  for  the  fact  that  atoms  spent 
only  about  10%  of  the  time  in  the  optical  potential  (the  rest  of  the  time  was  spent  in  the  pure 
magnetic  potential),  we  arrive  again  at  a  heating  rate  of  about  100  nK/s. 

Because  of  the  small  trapping  volume,  and  the  easy  availability  of  laser  light,  optically  trapped 
condensates  can  be  easily  compressed  to  extremely  high  densities:  in  some  of  our  work,  condensate 
densities  as  high  as  3  x  10*'^cm~^  were  observed.  At  such  densities,  the  three-body  recombina¬ 
tion  rate  is  extremely  high,  limiting  the  use  of  such  dense  samples.  At  the  low-density  end,  we 
have  produced  condensates  with  several  million  atoms  at  densities  of  about  1  x  10*^  cm^^  in  a 
decompressed  magnetic  trap  with  a  mean  frequency  u;  =  2t  x  7  IIz.  This  wide  range  of  possible 
condensate  densities  is  impressive,  especially  considering  that  a  few  years  ago  there  were  many 
discussions  whether  the  “BEC  window*'  even  exists. 

Besides  the  new  .scientific  po.ssibiIities  discussed  in  the  following  paragraphs,  the  optical  trap 
may  become  an  important  tool  to  manipulate  and  transport  condensates.  It  may  serve  as  an 
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“optical  tweezers”  to  move  condensates  for  example  into  optical  or  microwave  cavities  or  close  to 
surfaces. 


3.  REVERSIBLE  FORMATION  OF  A  BOSE-EINSTEIN  CONDENSATE 

In  an  ordinary  cryostat,  the  experimenter  can  raise  and  lower  the  temperature  of  the  sample 
reversibly.  In  contrast,  evaporative  cooling  is  irreversible  due  to  the  loss  of  the  evaporated  atoms. 
Even  if  the  temperature  is  raised  again  by  (internal  or  external)  heating  of  the  sample,  the  number 
of  atoms  has  been  reduced  during  the  cooling  stage.  This  reflects  the  fact  that  the  trapped  atoms 
are  both  the  sample  and  the  ‘Svorking  fluid^  of  the  refrigerator. 

If  the  strength  of  the  trapping  potential  is  changed  adiabatically,  one  can  reach  higher  or  lower 
temperatures  reversibly  -  but  the  density  changes  simultaneously  in  such  a  way  that  the  pha.se 
space  density  is  invariant.  Thus  adiabatic  changes  in  the  strength  of  the  trapping  potential  cannot 
Y>e  used  to  cross  the  BEC  transition  line. 

However,  let  us  consider  that  one  could  contain  a  small  portion  of  a  trapped  atomic  sample 
in  a  vessel  which  is  kept  in  thermal  contact  with  the  remaining  atoms,  and  then  use  this  vessel  to 
compre.ss  the  small  sample.  This  would  raise  the  density  of  a  small  portion  of  the  trapped  gas. 
but  because  of  the  thermal  contact  with  a  larger  reservoir  of  atoms,  the  temperature  rise  would  be 
slight.  Furthermore,  if  such  a  compression  were  performed  slowly  (quasi-statically).  the  entropy  of 
the  gas  would  be  unchanged,  and  one  could  reversibly  (by  compressing  and  expanding  the  small 
vessel)  change  the  maximum  phase-space  density  of  the  gas.  .4  similar  effect  is  achieved  without 
such  a  vessel  by  slowly  changing  the  shape  of  the  trapping  potential.  This  process  conserves  entropy, 
but  can  change  the  local  pha.se  space  density. 

This  adiabatic  (i.e.  isentropic)  increase  in  phase-space  density  can  only  occur  in  the  presence  of 
collisions.  Without  collisional  thermalization.  a  slow  adiabatic  process  would  preserve  the  number 
of  atoms  in  each  quantum  state  (i.e.  the  Ehrenfest  notion  of  adiabaticity)  and  therefore  would  not 
change  the  phase  space  density.  Thus,  the  maximal  increase  in  phase-space  density  from  changes 
in  the  trapping  potential  will  occur  if  such  changes  are  made  slowly  with  respect  to  the  collisional 
(and  motional)  equilibration  times  of  the  trapped  gas. 

The  possibility  of  adiabatic  changes  of  the  ground  state  population  was  recently  shown  by 
Pinkse  and  collaborators.^^  For  the  power-law  potentials  considered  in  their  paper,  the  maximum 
increase  of  phase-space  density  is  limited  to  a  factor  of  20.  A  more  general  deformation  of  the  trap¬ 
ping  potential  can  increase  the  phase-space  density  by  an  arbitrary  factor,  and  an  implementation 
of  this  scheme  was  realized  using  a  combination  of  magnetic  and  optical  forces. Atoms  in  the 
magnetic  trap  acted  as  the  reservoir  while  a  small  fraction  of  the  cloud  was  compressed  to  higher 
densities  by  adding  a  focused  red-detuned  laser  beam  which  creates  a  narrow  potential  well  (a  “dim¬ 
ple”)  in  the  center  of  the  magnetic  trap.  Such  small  “kinks”  of  the  trapping  potential  cannot  be 
realized  with  ordinary  magnetic  traps  which  are  produced  by  relatively  large  and  distant  magnets. 

The  increase  of  the  phase  space  density  was  measured  by  cooling  clouds  to  various  temperatures 
slightly  above  the  phase  transition  temperature  Tc  and  determining  the  infrared  laser  power  required 
to  produce  a  small  Bose  condensate.^^  To  attain  the  highest  possible  increa.se  in  phase-space 
density,  the  infrared  light  was  increased  slowly.  Fig.  1  shows  the  increases  in  pha.se-space  density. 
A  maximum  increa.se  by  a  factor  of  50  was  obtained.  Condensates  were  observed  in  clouds  with 
temperatures  as  high  as  5  //K.  Further  increa.ses  were  hindered  by  limitations  in  laser  power  and 
by  limits  to  the  ramp-up  time  set  by  the  various  heating  and  loss  processes  in  the  deformed  trap. 

Further  deformation  of  the  potential  takes  the  sample  across  the  BEC  phase  transition.  The 
observed  condensate  fractions  were  considerably  smaller  than  the  predictions  of  a  model  assuming 
non-interacting  atoms,  demonstrating  the  importance  of  interaction  effects  on  the  formation  of  the 
condensate.  The  importance  of  interactions  is  usually  parametrized  by  a  dimensionless  parameter 

which  is  defined  as  the  ratio  of  the  mean  field  energy  at  zero  temperature  to  the  transition 
temperature  Tc.  For  a  harmonic  trapping  potential  with  mean  oscillation  frequency  a;,  the  mean 
field  energy  is  proportional  to  and  Tc  oc  u;,  leading  to  ?/  oc  Due  to  this  weak  dependence 
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Fig.  1.  Phase-space  density  increase  to  reach 
BEC  vs.  normalized  well-depth  (from  Ref.  20). 
Various  trap  settings  were  used  for  the  the  mean 
trapping  frequency  u;  and  the  beam  waist  wq’  = 
X  100  IIz,  wq  =  9 /mi  (triangles);  u;  =  2:r  x  100 
Hz.  u'o  =  18/mi  (open  circles);  and  uj  =  2:r  x  33 
IIz.  u’o  =  iSfim  (closed  circles).  The  solid  line 
gives  theoretical  prediction. 


(and  also  a  weak  dependence  on  the  number  of  trapped  atoms),  all  experiments  done  so  far  realized 
i]  parameters  in  the  range  of  0.3  to  0.4^^  and  did  not  show  evidence  for  interaction  effects. 

In  contrast,  the  ‘'dimple“  trap  is  characterized  by  two  trapping  frequencies,  one  for  the  broad 
magnetic  trapping  potential  which  determines  Tc,  and  one  for  the  bottom  of  the  optical  trapping 
potential  which  confines  the  condensate  and  determines  its  mean  field  energy.  Thus,  the  parameter 
1]  can  be  greatly  increased,  exceeding  the  values  reached  in  purely  magnetic  traps  by  a  factor  of 
about  hundred.  Fig.  2  shows  the  ballistic  expansion  of  a  condensate  which  is  clearly  faster  than  the 
expansion  of  the  thermal  cloud,  which  means  that,  even  for  small  condensate  fractions,  the  mean 
field  energy  can  exceed  K'bT- 


Fig.  2.  In  the  dimple  trap  interaction  effects  are 
greatly  enhanced  and  the  meanfield  energy  of  the 
condensate  can  exceed  the  thermal  energy.  Shown 
is  the  ballistic  expansion  of  the  cloud  50  ms  after 
the  release  from  the  *Mimple”  trap.  The  disc¬ 
like  expansion  of  the  condensate  is  faster  than  the 
Isotropic  thermal  expansion. 

The  reversibility  of  crossing  the  BEC  phase  transition  was  demonstrated  by  preparing  a  mag¬ 
netically  trapped  cloud  just  above  Tc-  We  then  sinusoidally  modulated  the  power  of  the  infrared 
light  at  1  IIz.  between  zero  and  7  mW.  This  modulation  frequency  was  significantly  smaller  than 
the  magnetic  trap  frequencies  (^v  =  2;r  x  481  Iz  and  u;;  =  2;r  x  16  IIz).  These  low  frequencies  and 
a  large  optical  focus  (wq  =  18  //m)  were  used  to  minimize  trap  loss  due  to  inelastic  collisions. 

During  the  first  seven  condensation  cycles  the  condensate  fraction  oscillated  between  zero  and 
6  %.  as  shown  in  Fig.  3:  later  probing  showed  repeated  condensation  for  at  lea.st  15  cycles.  The 
peak  of  these  oscillations  decreased  slowly  in  time  due  to  a  small  heating  rate  resulting  from  slight 
instabilities  in  our  optical  setup. 

This  method  of  creating  condensates  provides  insight  into  their  formation,  which  was  recently 
studied  experimentally^^  and  theoretically.^^  For  example,  in  the  experiment  described  above,  the 
condensate  fraction  was  found  to  lag  about  70  ms  behind  the  modulation  of  the  la.ser  power,  which 
is  a  measure  for  the  formation  time.  Note  that  this  formation  time  is  much  smaller  than  the  1 
s  timescale  for  changing  the  trapping  potential,  thus  confirming  that  the  condensation  was  done 
quasi-statically. 

Compared  to  evaporative  cooling  in  magnetic  traps,  deformation  of  the  trapping  potential  offers 
several  advantages  for  the  study  of  the  BEC  pha.se  transition.  For  example,  in  a  previous  study 
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Fig.  3.  Adiabatic  cycling  through  the  phase  tran¬ 
sition  (from  Ref.  20).  Shown  is  the  condensate 
fraction  (a)  and  the  temperature  (b)  versus  time 
for  the  ca,se  of  a  modulated  infrared  beam  (closed 
circles),  an  infrared  beam  ramped  up  to  a  con¬ 
stant  power  (open  triangles),  and  no  infrared  light 
(open  squares).  The  solid  lines  are  guides  to  the 
eye. 


of  the  BEC  formation  process,  a  sudden  spatial  truncation  of  the  thermal  cloud  by  an  rf-sweep 
was  used  to  quench  the  system  into  BEC.^^  To  clearly  interpret  the  results  of  such  a  truncation, 
one  either  has  to  model  the  process  of  rf-truncation,  or  use  a  very  low  density  cloud  to  ensure  a 
separation  of  the  time  scales  for  truncation  and  thermalization.  In  contrast,  the  deformation  of  the 
trap  using  infrared  light  can  be  applied  in  microseconds  and  its  characterization  is  straightforward. 
In  one  such  experiment,  we  switched  on  the  infrared  light  instantly,  and  observed  condensation 
on  time  scales  much  faster  than  the  oscillation  periods  in  the  magnetic  trap  and  along  the  weakly 
confining  axis  of  the  optical  trap,  as  shown  in  Fig.  4.  The  resulting  condensates  showed  striations 
in  time-of“flight  images,  indicating  that  the  condensates  formed  into  excited  states  of  the  deformed 
potential-  Such  studies  of  shock-condensation  might  give  new  insight  into  the  formation  of  quasi¬ 
condensates  and  condensation  into  excited  states.^ 


5 

10  ; 

20 

30 

50  >  r 

70 

Fig.  4.  Formation  of  the  condensate  after  a  sudden  switch  on  of  the  IR  la.ser.  Shown  is  the  cloud 
after  a  ballistic  expansion  of  50  ms  for  various  dwell  times  in  the  ‘"dimple”  trap  (numbers  given  in 
ms).  A  small  condensate  fraction  is  already  visible  after  10  ms.  During  the  first  50  ms  of  formation 
the  condensate  shows  gross  striations  which  disappear  at  longer  dwell  times.  This  type  of  studies 
could  be  suitable  for  observing  the  formation  of  quasi-condensates,  where  phase  fluctuations  might 
lead  to  striations  in  the  ballistic  expansion. 


Also,  deformation  of  the  trapping  potential  greatly  extended  the  range  of  transition  temper¬ 
atures  which  we  can  easily  access.  The  compression  of  the  thermal  cloud  by  the  light  potential 
resulted  in  extreme  critical  densities  Uc  ^  b  x  10^’^ cm”^  and  critical  temperatures  Tc  ^  5  fiK,  At 
low  densities,  crossing  the  BEC  transition  in  a  purely  magnetic  trap  is  hindered  by  ineffectiveness 
of  rf-evaporation  (caused  by  the  sagging  of  the  cloud  due  to  gravity)  and  the  technical  difficulties 
of  controlling  very  small  trap  depth  with  rf-  and  magnetic  fields.  However,  using  our  method  of 
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deforming  the  potential,  we  could  shine  the  infrared  laser  into  a  decompressed  magnetic  trap  =  7 
IIz)  and  cross  the  transition  at  ^  2  x  10^^ cm*”’^  and  Tc  ^  100  nK. 

Reversible  cooling  by  adiabatic  deformation  of  the  trapping  potential  is  achieved  by  a  modi¬ 
fication  of  the  density  of  states.  When  a  potential  dimple  is  added  to  a  broad  harmonic  oscillator 
potential,  it  lowers  the  energy  of  the  ground  state,  but  affects  very  few  of  the  excited  states  of  the 
system.  As  the  well-depth  is  increased,  the  ground  state  energy  reaches  the  chemical  potential,  and 
DEC  occurs.  An  alternative  way  of  describing  this  process  is  to  first  compress  the  whole  cloud  by 
increasing  the  harmonic  oscillator  frequency  to  the  value  at  the  bottom  of  the  dimple,  and  then 
opening  up  the  potential  to  the  original  broad  harmonic  oscillator  potential  outside  a  small  central 
area.  The  last  step  increa.ses  the  density  of  states  at  high  energy,  and  as  a  result,  elastic  collisions 
preferably  populate  these  .states.  This  is  very  similar  to  evaporative  cooling  where  rf-induced  spin- 
flips  couple  higher  lying  states  to  the  continuum  leading  to  a  net  collisional  transfer  of  atoms  to 
these  states.  The  major  difference  is  that  in  evaporative  cooling  no  equilibrium  is  established  be¬ 
tween  the  evaporated  atoms  and  the  trapped  atoms,  whereas  in  the  ca.se  of  the  deformed  potential, 
the  whole  cloud  stays  in  thermal  equilibrium. 

4.  OBSERVATION  OF  FESHBACH  RESONANCES  IN  A  BOSE-EINSTEIN 

CONDENSATE 

All  the  essential  properties  of  Bose  condensed  systems  -  the  formation  and  shape  of  the  conden- 
.sate.  the  nature  of  its  collective  excitations  and  statistical  fluctuations,  the  formation  and  dynamics 
of  solitons  and  vortices  -  are  determined  by  the  strength  of  the  atomic  interactions.  In  contra.st  to 
the  situation  for  superfiuid  helium,  the.se  interactions  are  weak,  permitting  the  phenomena  to  be 
theoretically  described  from  "first  principles**.  Furthermore,  in  atomic  gases  the  interactions  can  be 
altered,  for  instance  by  employing  different  species,  changing  the  atomic  density,  or.  as  we  discuss 
here,  merely  by  varying  a  magnetic  field. 

At  low  temperatures,  the  interaction  energy  in  a  cloud  is  proportional  to  the  density  and  a 
single  atomic  parameter,  the  scattering  length  a.  which  depends  on  the  quantum-mechanical  pha.se 
shift  in  an  elastic  collision.  It  has  been  predicted  that  the  scattering  length  can  be  modified  by 
applying  external  magnet optical'^^'^^  radio-frequency  (rf).‘^‘^  and  electric  fields.'^*  Those 
modifications  are  only  pronounced  near  zero-energy  resonances,  either  shape  resonances’^*  or  "Fes- 
hbach  re.sonances*‘.’^"  A  Feshbach  resonance  occurs,  when  a  qua.sibound  molecular  state  has  nearly 
zero  energy  and  couples  resonantly  to  the  free  state  of  the  colliding  atoms.  In  a  time-dependent 
picture,  the  two  atoms  are  transferred  to  the  quasibound  state,  "stick**  together  and  then  return 
to  an  unbound  state.  Such  a  resonance  strongly  affects  not  only  the  scattering  length  (elastic 
channel),  but  also  inelastic  processes  such  as  dipolar  relaxation^^*^'  and  three-body  recombination. 
Feshbach  resonances  have  so  far  been  studied  at  much  higher  energies’^^  by  varying  the  collision 
energy.  In  our  recent  work^'  we  have  shown  that  they  can  be  tuned  to  zero  energy  to  be  resonant 
for  ultracold  atoms.  The  different  magnetic  moments  of  the  free  and  qua.sibound  states  allowed  us 
to  tune  these  resonances  with  magnetic  fields,  and  as  a  result,  minute  changes  in  the  magnetic  field 
strongly  affected  the  properties  of  a  macroscopic  sy.stem. 

.Above  and  below  a  Feshbach  resonance,  the  scattering  length  a  covers  the  full  continuum  of 
positive  and  negative  values.  It  was  predicted  to  vary  dispersively  as  a  function  of  magnetic  field 


where  A  parameterizes  the  width  of  the  resonance  at  B  ^  Dq,  and  a  is  the  scattering  length  away 
from  the  re.sonance.  Feshbach  re.sonances  should  allow  the  realization  of  condensates  over  a  wide 
range  of  interaction  strengths.  By  setting  a  ^  0,  one  can  create  a  condensate  with  essentially 
non-interacting  atoms,  and  by  setting  a  <  0  one  can  make  the  system  unstable  and  ob.serve  its 
collapse.  Rapid  tuning  of  an  external  magnetic  field  around  a  Feshbach  re.sonance  will  lead  to 
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sudden  changes  of  the  scattering  length.  This  opens  the  door  to  studies  of  new  dynamical  effects 
such  as  novel  forms  of  collective  oscillations  or  the  sudden  collapse  of  a  large  condensate  when 
the  scattering  length  is  switched  from  positive  to  negative.^^  Feshbach  resoncances  have  recently 
been  also  observed  in  rubidium,  in  photoassociation  spectroscopy  of  laser-cooled  atoms.^^  and  as  a 
resonance  in  the  elastic  scattering  cross  section  of  magnetically  trapped  atoms.*^ 

Theoretical  calculations  predicted  Feshbach  resonances  for  sodium  only  for  hyperfine  states 
which  cannot  be  magnetically  trapped. The  optical  trap  described  above  was  therefore  in- 
dispensible  for  this  study.  After  optically  trapped  Bose-Einstein  condensates  were  produced  as 
described  in  Section  2.  they  could  be  transfered  to  the  mf  =  H-l  state  with  an  adiabatic  rf~sweep. 
and  were  then  exposed  to  a  bias  field  up  to  ~  1200  G.  The  in  homogeneity  of  the  bias  field  destabi¬ 
lized  the  optical  trap  in  the  axial  direction  and  therefore  required  the  use  of  twoadditonal  "‘endcaps 
of  far-off- resonant  blue-detuned  laser  light.^' 

When  the  magnetic  field  is  swept  across  a  Feshbach  resonance  one  would  expect  to  lose  a 
condensate  due  to  an  enhanced  rate  of  inelastic  collisions  (caused  either  by  the  collapse  in  the 
region  of  negative  scattering  length  or  by  an  enhanced  rate  coefficient  for  inelastic  collisions).  This 
was  used  to  locate  three  resonances  near  853  G.  907  G.  and  1195  G  by  the  sudden  onset  of  trap 
loss. 

In  order  to  measure  the  variation  of  the  scattering  length  around  the  resonance  near  90 r  G.  we 
determined  the  interaction  energy  of  a  trapped  condensate.  This  was  done  by  suddenly  switching  off 
the  trap,  allowing  the  stored  interaction  energy  to  be  converted  into  the  kinetic  energy  of  a  freely 
expanding  condensate  and  measuring  it  in  time-of-flight  absorption  imaging.  The  interaction 
energy  is  proportional  to  the  scattering  length  and  the  average  density  of  the  condensate.  Both 
quantities  can  be  directly  extracted  from  time-of-flight  images,  and  the  results  for  the  scattering 
length  are  shown  in  Fig.  5  together  with  the  theoretical  prediction*^  of  a  A  =  1  G  wide  resonance. 
The  data  clearly  displays  the  predicted  dispersive  shape  and  shows  evidence  for  a  variation  in  the 
scattering  length  by  more  than  a  factor  of  ten. 


Magnetic  Field  (G) 


Fig.  5.  Observation  of  the  Feshbach  resonance 
at  907  G  using  time-of-flight  absorption  imag¬ 
ing  (from  Ref.  37).  Fig.  a)  shows  the  num¬ 
ber  of  atoms  in  the  condensate  versus  magnetic 
field.  Field  values  above  the  resonance  were 
reached  by  quickly  crossing  the  resonance  from 
below  and  then  approaching  from  above  slowly. 
Fig.  b)  shows  the  normalized  scattering  length 
a/ a  a  calculated  from  the  relea.sed  en¬ 

ergy,  together  with  the  predicted  shape  (Eq.  1, 
solid  line).  The  values  of  the  magnetic  field  in 
the  upper  scan  relative  to  the  lower  one  have  an 
uncertainty  of  less  than  0.5  G. 


The  trap  losses  observed  around  the  Feshbach  resonances  merit  further  study  since  they  might 
impose  practical  limits  on  the  possibilities  for  varying  the  scattering  length.  An  increase  of  the 
dipolar  relaxation  rate  near  Feshbach  resonances  has  been  predicted, but  for  atoms  in  the 
lowest  hyperfine  .state  no  such  inelastic  binary  collisions  are  possible.  Therefore,  the  observed  trap 
loss  for  the  resonances  in  the  \F  -  1,  inp  -  H-l)  state  is  probably  due  to  three-body  collisions.  In 


J.  St  eager  et  aL 

this  case  the  loss  rate  is  characterized  by  the  coefficient  A'3,  defined  as  N /N  =  -A’3  (n^).  So  far. 
there  is  no  theoretical  work  on  A3  near  a  Feshbach  resonance.  An  analysis  based  on  Fig.  5  shows 
that  A'3  increased  on  both  sides  of  the  resonance,  because  the  loss  rate  increased  while  the  density 
decreased  or  stayed  constant.  The  inferred  increase  of  A'3  is  fairly  dramatic  (up  to  five  orders  of 
magnitude).  The  trap  loss  prevented  the  study  of  the  region  with  negative  scattering  length:  the 
condensate  was  lost  before  this  region  could  be  reached.  A  lifetime  of  the  condensate  of  at  least  one 
oscillation  period  in  the  trap  is  required  to  observe  the  collapse  of  an  unstable  condensate.  For  the 
sodium  resonances  presented  here,  this  requires  lower  density  condensates  and  therefore  an  optical 
trap  with  a  larger  volume.  Whether  the  observed  enhancement  of  the  three-body  recombination 
rate  is  generic  for  any  Feshbach  resonance  or  only  specific  to  sodium  remains  to  be  seen  -  both 
experimental  and  theoretical  work  along  these  lines  is  in  progress. 

Our  observation  of  the  dispersive  variation  of  the  scattering  length  confirms  the  theoretical 
predictions  about  “tunability"  of  the  scattering  length  with  the  prospect  of  ‘designing*'  atomic 
quantum  gases  with  novel  properties:  e.g.  to  create  ideal  Bose-Einstein  condensates  with  nearly 
zero  .scattering  length  and  to  obtain  a  detailed  picture  of  the  collapse  of  a  condensate  with  negative 
scattering  length  which  is  so  far  not  fully  understood.  Tuning  the  scattering  length  can  also  be 
used  to  vary  interactions  between  different  specie.s'^^  and  thus  control  the  pha.se  diagram  of  multi- 
component  condensates,  possibly  switching  from  interpenetrating  superfluids  to  phase  separation.’*’^ 
Feshbach  resonances  may  also  be  important  in  atom  optics,  for  modifying  the  atomic  interactions  in 
an  atom  laser,  or  more  generally,  for  controlling  non-linear  coefficients  in  atom  optics  with  coherent 
beams  of  atoms. 


5.  SPINOR  BOSE-EINSTEIN  CONDENSATES  IN  OPTICAL  TRAPS 

Bose-Einstein  condensates  in  magnetic  traps  are  described  by  a  scalar  order  parameter  sim¬ 
ilar  to  the  spinle.ss  superfluid  **IIe.  Even  though  alkali  atoms  have  angular  momentum,  the  spin 
orientation  is  not  a  degree  of  freedom  because  the  atomic  spin  adiabatically  follows  the  magnetic 
field.  Spin  flips  lead  to  untrapped  states  and  are  therefore  a  loss  proce.ss.  In  contra.st.  an  optical 
trap  confines  atoms  in  all  spin  states,  opening  the  possibility  to  .study  spinor  condensates  which 
have  the  spin  orientation  as  a  new  degree  of  freedom.  They  represent  a  system  with  a  vector  order 
parameter  instead  of  a  scalar.  A  variety  of  new  phenomena  is  predicted,  *  *”  *^  such  as  spin  textures, 
propagation  of  spin  waves  and  coupling  between  superfluid  flow  and  atomic  spin. 

Other  systems  which  go  beyond  the  description  with  a  single  scalar  order  parameter  are  con- 
den.sates  of  two  different  hyperfine  states  of  ^'  Rb  confined  in  magnetic  traps.  Recent  experimental 
.studies  have  explored  the  spatial  separation  of  the  two  components'*^' *^  and  their  relative  phase.  *^ 
Several  theoretical  papers  de.scribe  their  structure*'^* and  their  collective  excitations."'*^''*^ 

Compared  to  these  two-component  condensates,  spinor  condensates  have  .several  new  features 
including  the  multi-component  vector  character  of  the  order  parameter  and  the  changed  role  of  spin 
relaxation  collisions  which  allow  for  population  exchange  among  hyperfine  states  without  trap  loss. 
In  contrast,  in  the  ^*Rb  experiments  trap  loss  due  to  spin  relaxation  severely  limits  the  lifetime. 

We  have  recently  studied  the  equilibrium  .state  of  spinor  condensates  in  an  optical  trap.^^  In 
an  F=1  spinor  condensate  subject  to  spin  relaxation,  two  m/r  =  0  atoms  can  collide  and  produce  a 
mjr  =  -^i  and  a  m/7  =  — 1  atom  and  vice  versa.  The  mo.st  dramatic  effect  is  seen  when  we  start  out 
with  a  condensate  in  a  pure  m=0  state.  Depending  on  the  external  magnetic  field,  we  observed  the 
formation  of  three  domains  of  m/r=  and  —1  atoms. 

The  experimental  study  of  spinor  condensates  requires  techniques  to  selectively  prepare  and 
probe  condensates  in  arbitrary  hyperfine  states.  The  first  step  was  the  preparation  of  n?/r  =  -1 
condensates  in  an  optical  trap  as  de.scribed  above.  Arbitrary  populations  of  the  three  hyperfine 
states  were  prepared  using  rf-transitions  (Landau-Zener  sweeps)  similar  to  our  earlier  work.'*  At 
a  bias  field  of  about  40  G  the  transitions  from  7?2/r  =  -l  to  m/r  =  0  and  from  mjr  =  0  to  77?/r  =  -fl 
differ  in  frequency  by  about  O.OMIIz  due  to  the  quadratic  Zeeman  shift  and  they  could  be  driven 
separately.  Without  this  frequency  difference  rf-fields  would  only  rotate  the  spin  vector  and  it 
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would  not  be  possible  to  produce  pure  mp  —  0  condensates.  After  the  spin  preparation,  a  bias  field 
Dq  and  a  field  gradient  D'  were  applied  for  a  variable  amount  of  time  (as  long  as  30  s).  during 
which  the  atoms  relaxed  towards  their  equilibrium  distribution.  The  bias  field  was  varied  between 
20mG  and  oOOniG  and  the  gradients  between  ±70mG/cm. 

Both  the  spatial  and  hyperfine  distributions  were  observed  by  time-of-flight  absorption  imag¬ 
ing.  For  this,  the  optical  trap  was  suddenly  switched  off  and  the  atoms  were  allowed  to  expand. 
Due  to  the  large  aspect  ratio,  the  expansion  was  almost  purely  in  the  radial  directions.  A  magnetic 
field  gradient,  which  was  applied  after  5  ms  time-of-ffight.  acted  as  a  Stern-Gerlach  filter  and  sep¬ 
arated  the  three  spin  components.  In  this  manner,  a  single  time-of-fiight  image  provided  both  a 
spatial  and  a  spin-state  description  of  the  trapped  cloud.  After  a  total  time-of-flight  of  25  ms  the 
atoms  were  pumped  into  the  F  =  2  hyperfine  state  and  observed  using  the  mjr=:-\-2  to  m/r  =  -l-.3 
cycling  transition.  This  technique  assured  the  same  transition  strength  for  atoms  originating  from 
different  spin  states. 
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Fig.  6.  Absorption  images  of  spinor  condensates  released  from  the  optical  trap,  after  25  msec 
of  ballistic  expansion  (from  Ref.  60).  The  different  mjr  states  were  separated  by  an  axial  field 
gradient.  The  images  taken  after  various  dwell  times  in  the  trap  show  how  the  atoms  evolved  into 
the  same  equilibrium  distribution  although  they  were  initially  prepared  in  a  pure  mr  =  0  .state 
(upper  row)  or  in  equally  populated  7iif  =  il  states  (lower  row).  The  bias  field  during  the  dwell 
time  was  Dq  =  20  mG  and  the  gradient  was  Z?'  =  11  niG/cm.  The  height  of  the  images  is  2.7  mm. 

Fig.  6  shows  a  sequence  of  images  with  different  dwell  times  in  the  optical  trap.  Starting  with 
the  pure  Wf  =  0  component  (upper  series)  or  with  a  50-50  mixture  of  the  inf  =  ±i  components 
(lower  series),  the  same  equilibrium  distribution  was  reached. 

Comparison  of  images  like  those  in  Fig.  6  to  a  theoretical  model  revealed  that  the  spin- 
dependent  interaction  cFi  •  F2  between  two  sodium  atoms  in  the  F  =  1  state  is  anti-ferromagnetic 
(i.e.  c  >  0),  Furthermore,  the  experimental  results  showed  clear  evidence  for  the  miscibility  of 
Wf  =  ~1  and  mf  =  -fl  components  and  the  immiscibility  of  inf  =  dbl  and  lUf  =  0.  This  opens 
the  possibility'  for  detailed  studies  of  miscible  and  immiscible  multi-component  condensates. 

The  anti-ferromagnetic  energy  was  estimated  to  be  2.5  iiK  in  our  condensates.  It  is  interesting 
to  note  that  the  optically  trapped  samples  in  which  the  domains  were  observed  were  at  a  temper¬ 
ature  of  the  order  of  100  nK.  far  larger  than  the  anti-ferromagnetic  energy.  Thus  the  formation  of 
spin  domains  occurs  only  in  a  Bose-Einstein  condensate.  This  illustrates  a  new  regime  acce.ssible 
with  Bose-Einstein  condensates:  a  few  nK  energy  are  responsible  for  spin  domains  of  a  gaseous 
Bose-condensed  cloud. 

The  dynamics  of  domain  formation  deserves  further  study.  For  example  we  have  placed  a 
two-component  condensate  in  a  state  with  higher  energy  than  the  groundstate.  representing  a 
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non-equilibriuni  configuration  which  is  nietastable  for  many  seconds.  We  are  currently  studying  the 
conditions  for  this  inetastability  and  the  mechanisms  by  which  the  equilibrium  finally  is  approached. 

6.  THE  SHAPE  OF  MAGNETICALLY  TRAPPED  BOSE-EINSTEIN 

CONDENSATES 

Finally,  we  would  like  to  present  some  unpublished  results  from  our  work  with  magnetically 
trapped  condensates:  in  situ  density  profiles  of  atomic  clouds  above  and  below  the  BEC  phase 
transition  which  illustrate  the  onset  of  a  spatial  condensation  of  atoms,  and  measurements  on 
condensates  released  from  the  magnetic  trap  and  allowed  to  freely  expand,  which  confirm  the 
mean^field  description  of  such  condensates. 

Bose-Einstein  condensation  in  a  homogeneous  system  is  appropriately  described  as  "conden¬ 
sation  in  momentum  space''^^  because  it  does  not  lead  to  spatial  separation  between  the  normal 
component  and  the  condensate.  However,  in  any  inhomogeneous  potential  -  i.e.  atom  traps  -  the 
condensate  and  the  normal  fraction  are  spatially  separated.^^’^^  This  occurs  because  the  condensate 
atoms  settle  into  the  areas  of  lowest  potential  energy  where  the  ground  state  lies.  In  situ  images  of 
such  a  separation  were  first  taken  with  alkali  atom  gases  using  the  dark-ground  imaging  technique, 
which  clearly  distinguished  the  high-density  of  the  condensate  from  the  background  thermal  cloud 
in  a  magnetically  trapped  sample.^^  Since  then,  the  phase-contrast  technique  has  been  used  for  in 
situ  studies  of  the  condensate,  using  the  non-destructive  character  of  such  non-resonant  imaging 
to  image  the  condensate  repeatedly  and  measure  dynamical  properties  such  as  the  propagation 
of  Bogoliubov  sound.^^  the  formation  of  a  condensate.^^  and  collective  excitations  at  non-zero 
temperatures.'^'*  In  such  images,  the  signal  which  is  obtained  is  linear,  rather  than  quadratic,  in 
the  column  density  of  the  gaseous  sample,  and  thus  clear  images  with  a  large  dynamic  range  are 
obtained. 
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Fig.  7.  Bose-Einstein  condensation  in  coordi¬ 
nate  space.  Traces  show  axial  profiles  of  column 
densities  of  Bose  gases  from  phase-contrast  im¬ 
ages.  Trapping  frequencies  are  17(1)  IIz  axially 
and  230(20)  IIz  radially.  Temperatures  are  deter¬ 
mined  from  fits  to  the  thermal  tails.  As  the  tem¬ 
perature  is  lowered,  the  thermal  cloud  is  depleted, 
and  atoms  accumulate  in  the  Bose  condensate. 


In  Fig.  7,  we  present  a  series  of  profiles  of  column  densities  obtained  from  phase-contra.st 
images  of  clouds  above  and  below  the  BEC  phase  transition.  The  profiles  are  taken  along  the 
axis  in  which  the  condensate  is  longest.  As  in  our  other  work,  such  images  are  taken  with  light 
which  is  1.7  GIIz  red-detuned  from  the  D2  resonance  of  a  magnetically  trapped  sample  of  .sodium 
atoms  in  the  \F  ~  1,  m/r  =  —1)  h.yperfine  state.  While  such  images  are  non-destructive,  the  five 
traces  shown  in  the  figure  are  taken  from  different  clouds,  each  of  which  was  brought  to  a  different 
temperature.  One  can  clearly  see  that  as  the  temperature  is  lowered  from  1.8  pK  to  1.2  and 
below,  a  high-density  peak  arises  in  the  center  of  the  cloud,  corre.sponding  to  the  accumulation 
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of  condensate  atoms  at  the  bottom  of  the  confining  potential  well.  As  the  temperature  is  lowered 
further,  to  below  0.5  f.tK.  the  "thermal  tails"  of  the  central  condensate  peak  disappear,  indicating 
that  almost  all  the  remaining  atoms  form  a  Bose-Einstein  condensate,  which  contains  on  the  order 
of  10^  atoms.  The  width  of  the  condensate  (almost  300  //m  at  low  temperatures  in  this  figure) 
is  determined  by  its  mean-field  energy:  at  the  edge  of  the  condensate,  the  trap  potential  energy 
and  the  mean-field  energy  (chemical  potential)  of  the  condensate  are  equal.  In  comparison,  the 
non-interacting  ground  state  in  a  harmonic  well  with  a  trap  frequency  of  17  IIz  (as  in  the  figure) 
has  a  half-width  of  5  ^m  -  a  clear  indication  of  the  importance  of  interactions  in  the  .static  and 
dynamical  properties  of  such  condensates,  and  a  clear  justification  for  the  use  of  the  Thomas-Fermi 
approximation. 

The  temperature  of  each  cloud,  as  indicated  on  each  trace,  is  determined  by  fitting  a  Gaussian 
function  to  the  wings  of  the  density  profile,  away  from  the  central  condensate  peak.  The  width  of 
this  Gaussian  can  be  combined  with  the  trapping  frequency  in  the  axial  direction  (which  can  be 
determined  precisely  by  in  situ  measurements)  to  obtain  the  temperature.  However,  as  indicated 
for  the  last  two  traces  in  the  figure,  temperatures  which  are  lower  than  the  chemical  potential  are 
difficult  to  measure  by  this  method,  because  the  extent  of  the  thermal  cloud  and  of  the  condensate 
are  similar,  and  furthermore  the  non-condensed  fraction  at  such  low  temperatures  is  extremely 
small.  This  poses  a  limitation  on  studying  low-temperature  properties  of  condensates  (as  discussed 
in  Section  2). 

This  problem  may  be  somewhat  obviated  by  measuring  temperatures  in  anisotropic  traps  by 
1 1 me-of- flight  techniques.  Since  condensates  from  such  traps  expand  primarily  in  only  one  (pancake 
condensates)  or  two  (cigar-shaped  condensates)  directions  while  the  thermal  cloud  expands  isotrop¬ 
ically.  a  clearer  spatial  separation  of  the  two  components  occurs  along  the  axis  of  slow  condensate 
expansion,  allowing  for  a  temperature  measurement. 

This  anisotropic  expansion  of  Bose-Einstein  condensates  from  harmonic  magnetic  traps  has 
been  investigated  experimentally  and  theoretically.^^  The  anisotropy  was  one  of  the  distinguishing 
characteristics  of  Bose-Einstein  condensates  which  allowed  for  their  identification  in  early  BEC  ex¬ 
periments.  and  allows  for  the  clear  separation  between  the  condensed  and  non-condensed  fractions. 
The  shape  of  the  condensate  in  time-of- flight  arises  from  the  anisotropy  of  a  trapped  condensate  for 
two  reasons.  First,  the  shape  of  the  trapped  condensate  implies  an  anisotropic  momentum  distribu¬ 
tion  in  the  trap.  thus,  in  a  interaction-free  time-of-flight.  this  momentum  distribution  should  yield 
an  anisotropic  expanding  cloud.  However,  in  the  presence  of  interactions,  this  kinetic  energy  effect 
is  much  smaller  than  the  mean-field  energy  of  the  condensate;  this  energy  is  released  preferentially 
along  the  narrowest  directions  of  the  trapped  condensate  and  is  converted  to  kinetic  energy.  For 
large  condensates,  this  expansion  is  well-described  by  the  Gross- Pitaevskii  equation  solved  in  the 
Thomas-Fermi  approximation^* ^  wherein  the  evolution  is  described  by  a  simple  scaling  solution 
of  the  original  parabolic  density  profile  in  the  harmonic  trap. 


Fig.  8.  Aspect  ratio  (axial  to  radial  widths)  of 
expanding  Bose  condensates.  Data  collected  for 
times  less  than  10  ms  were  taken  in  pha.se-contra.st 
imaging  due  to  the  high  optical  density.  Other 
data  were  taken  with  absorption  imaging.  The 
dashed  line  indicates  the  ratio  of  trapping  fre¬ 
quencies  of  248(1)  Hz  radially  and  16.23(3)  Hz 
axially.  The  solid  line  gives  the  prediction  of  Ref. 
69. 


We  probed  this  free-evolution  using  a  combination  of  phase-contrast  (short  time-of-flight) 


J.  Steiiger  et  al. 


and  absorption  (long  time-of-flight)  imaging.  The  measured  aspect  ratios  of  a  freely-expanding 
condensate  are  presented  in  Fig.  8.  The  data  clearly  show  two  steps  in  the  expansion.  First,  during 
the  first  millisecond  or  so.  the  atoms  are  accelerated  radially  by  the  stored  mean-field  energy.  Then, 
the  atoms  travel  ballistically,  expanding  the  condensate  in  the  radial  direction  at  a  constant  velocity, 
while  the  axial  expansion  is  nearly  absent  due  to  the  large  original  axial  extent  of  the  cloud.  The 
final  stage  of  the  expansion,  where  the  cloud  reaches  a  final,  asymptotic  aspect  ratio,  occurs  only 
at  times  larger  than  our  accessible  range  of  times-of-Wight.  The  data  are  presented  along  with  a 
predicted  evolution^^  for  which  the  only  input  is  the  axial  and  radial  trap  frequencies  which  we 
measured  precisely  with  repeated  in  s/tn  imaging.  The  agreement  between  theory  and  experiment  is 
excellent,  and  provides  convincing  proof  for  the  validity  of  the  Gross-Pitaevskii  equation  to  describe 
large-amplitude  dynamics. 
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Bose  Einstein  condensates  of  dilute  atomic  gases,  characterized  by  a  macroscopic  pop¬ 
ulation  of  the  quantum  mechanical  ground  state,  are  a  new,  weakly  interacting  quantum 
fluid^’  In  most  experiments  condensates  in  a  single  weak  field  seeking  state  are  magnet¬ 
ically  trapped.  These  condensates  can  be  described  by  a  scalar  order  parameter  similar 
to  the  spinless  superfluid  ^He.  Even  though  alkali  atoms  have  angular  niomentum,  the 
spin  orientation  is  not  a  degree  of  freedom  because  spin  flips  lead  to  un trapped  states 
and  are  therefore  a  loss  process.  In  contrast,  the  recently  realized  optical  trap  for  sodium 
condensates'^  confines  atoms  independently  of  their  spin  orientation.  This  opens  the  possi¬ 
bility  to  study  spinor  condensates  which  represent  a  system  with  a  vector  order  parameter 
instead  of  a  scalar.  Here  we  report  a  study  of  the  equilibrium  state  of  spinor  condensates  in 
an  optical  trap.  The  freedom  of  spin  orientation  leads  to  the  formation  of  spin  domains  in 
an  external  magnetic  field.  The  structure  of  these  domains  are  illustrated  in  spin  domain 
diagrams.  Combinations  of  both  miscible  and  immiscible  spin  components  were  realized. 

A  variety  of  new  phenomena  is  predicted^'  •  for  spinor  condensates,  such  as  spin  textures,  propa¬ 
gation  of  spin  waves  and  coupling  between  superfluid  flow  and  atomic  spin.  To  date  such  effects  could 
only  be  studied  in  superfluid  ^He,  which  can  be  described  by  Bose  Einstein  condensation  of  Cooper 
pairs  of  quasi  particles  having  both  spin  and  orbital  angular  momentum’*.  Compared  to  the  strongly 
interacting  ^He,  the  properties  of  weakly  interacting  Bose  Einstein  condensates  of  alkali  gases  can  be 
calculated  by  mean  field  theories  in  a  much  more  straightforward  and  simple  way. 

Other  systems  which  go  beyond  the  description  with  a  single  scalar  order  parameter  are  condensates 
of  two  different  hyperfine  states  of  ^'Rb  confined  in  magnetic  traps.  Recent  experimental  studies  have 
explored  the  spatial  separation  of  the  two  components^'’-  and  their  relative  phase^^  Several  theoretical 
papers  de.scribe  their  structure^"-  and  their  collective  excitations^'^- 

Compared  to  these  two  component  condensates,  spinor  condensates  have  several  new  features  in¬ 
cluding  the  vector  character  of  the  order  parameter  and  the  changed  role  of  spin  relaxation  collisions 
which  allow  for  population  exchange  among  hyperfine  states  without  trap  loss.  In  contrast,  for  ^‘Rb 
experiments  trap  loss  due  to  spin  relaxation  severely  limits  the  lifetime. 

We  consider  an  F=1  spinor  condensate  subject  to  spin  relaxation,  in  which  two  /777:’  =  0  atoms  can 
collide  and  produce  an  777/'  = +  1  and  an  nip  =  —\  atom  and  vice  versa.  We  investigate  the  distribution 
of  hyperfine  states  and  the  spatial  distribution  in  equilibrium  assuming  conservation  of  the  total  spin. 

The  ground  state  spinor  wave  function  is  found  by  minimizing  the  free  energy^ 

K  =  J  d^rn  +  ^  +  ^{Ff  +  E.-.  -  />o(Er>]  ,  (1) 

where  kinetic  energy  terms  are  neglected  in  the  Thomas- Fermi  approximation  which  is  valid  as  long 
as  the  dimension  of  spin  domains  (typically  50  fim)  is  larger  then  the  penetration  depth^*  (typically 
1  //777).  V  is  the  trapping  potential,  n  is  the  density,  F  is  the  angular  momentum  per  atom,  and 
is  the  Zeeman  energy  in  an  external  magnetic  field.  The  Lagrange  multiplier  po  accounts  for  the  total 
spin  conservation.  The  mean  field  energy  in  Eqn.  (1)  consists  of  a  spin  independent  part  proportional 
to  Co  and  a  spin  dependent  part  proportional  to  C2{F)‘.  The  coefficients  cq  and  co  are  related  to  the 
scattering  lengths  ao  and  ao  for  two  colliding  atoms  with  total  angular  momentum  Ftot  =  0  or  Ftot  =  2 
by  Co  =  47rTra/M  and  C2  =  47TtrAa/M  with  a  =  (2a2  +  ao)/3,  Aa  =  (02  —  uo)/3,  and  M  for  the 
atomic  mass'".  The  spin  dependent  interaction  originates  from  the  term  CoFi  *  F2  in  the  interaction  of 
two  atoms,  which  is  ferromagnetic  for  C2  <  0  and  anti  ferromagnetic  for  C2  >  0. 
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In  the  Bogoliiibov  approach  the  many  body  ground  state  wave  function  is  represented  by  the  spinor 
wave  function 

=  v/i^<('")  =  \/^  (C+(r),Co(0>C-(0),  (2) 

where  denote  the  amplitudes  for  the  +1,0,— 1  states,  respectively,  and  |4‘|‘  =  1. 

The  Zeeman  energy  E^e  is  given  by 

+  ^olCol’  +  =  Eo  -p{Ez)  +  g  (F;).  (3) 

E+.  Eq.  E^  are  the  Zeeman  energies  of  the  =  +1, 0,  -1  states,  2^  =  E+  +  E_  ~  2Eo  is  the  Zeeman 
energy  difference  in  a  spinflip  collision,  and  2p  =  The  Eq  term  can  be  included  in  the  trapping 

potential  V.  The  parameter  p  can  be  combined  with  the  Lagrange  multiplier  pQ  to  give  p  =  p  +  po- 

In  the  following  we  determine  the  spinor  which  minimizes  the  spin  dependent  part  Kg  of  the  free 
energy: 

K,=c{Ff-p{F,)  +  q{F:-),  (4) 

where  c  =  The  minimization  of  Eqn.  (4)  for  different  values  of  the  parameters  c,  p,  and  q  is 

straightforward,  and  is  shown  graphically  in  the  form  of  spin-domain  diagrams  in  Fig.  1, 

Experimentally,  the  values  of  c,  p.  and  q  can  be  varied  arbitrarily,  representing  any  region  of  the 
spin  domain  diagram.  The  magnitude  (but  not  the  sign)  of  the  coefficient  c  is  varied  by  changing  the 
density  n,  either  by  changing  the  trapping  potential,  or  by  studying  condensates  with  different  numbers 
of  atoms.  In  this  study,  the  axial  length  of  the  trapped  condensate  is  more  than  60  times  larger  than  its 
radial  size,  and  thus  we  consider  the  system  one  dimensional,  and  integrate  over  the  radial  coordinates, 
obtaining  n  =  2no/3  where  uq  is  the  density  at  the  radial  center.  This  integration  assumes  a  parabolic 
density  profile  within  the  Thomas  Fermi  approximation.  The  value  of  q  can  be  changed  by  applying  a 
weak  external  bias  field  Eq:  q  then  corresponds  to  the  quadratic  Zeeman  shift  q  =  qB^.  The  coefficient  p 
arises  both  from  the  linear  Zeeman  shift  and  from  the  Lagrange  multiplier  po  which  is  determined  by  the 
total  spin  of  the  system.  For  a  system  with  zero  total  spin  in  a  homogenous  bias  field  Eo,  po  cancels  the 
linear  Zeeman  shift  due  to  Eq,  yielding  p  =  0.  Positive  (negative)  values  of  p  are  achieved  for  condensates 
with  a  positive  (negative)  overall  spin.  Finally,  the  coefficients  can  be  made  to  vary  spatially  across  the 
condensate.  In  particular,  applying  a  field  gradient  E'  along  the  axis  of  the  trapped  condensate  causes 
p  to  vary  along  the  conden.sate  length.  For  a  condensate  with  zero  total  spin,  p  =  pE'c  where  c  is  the 
axial  coordinate  with  c  =  0  at  the  center  of  the  condensate.  Thus,  the  condensate  samples  a  vertical 
line  in  the  spin  domain  diagrams  of  Fig.  1.  The  center  of  this  line  lies  at  p  =:  0,  and  its  length  is  given 
by  the  condensate  length  scaled  by  pE'. 

The  experimental  study  of  spinor  condensates  required  techniques  to  selectively  prepare  and  probe 
condensates  in  arbitrary  hyperfine  states.  Spinor  condensates  were  prepared  in  several  steps.  Lciser 
cooling  and  evaporative  cooling  were  used  to  produce  sodium  condensates  in  the  z??/’ =  —1  state  in  a 
cloverleaf  magnetic  trap“^.  The  condensates  were  then  tran.sfered  into  an  optical  dipole  trap  consisting  of 
a  single  focused  infrared  laser  bearn^  Arbitrary  populations  of  the  three  hyperfine  states  were  prepared 
using  rf  transitions^.  After  the  spin  preparation,  a  bias  field  Eq  and  a  field  gradient  E'  were  applied  for 
a  variable  amount  of  time  (as  long  as  30  s),  during  which  the  atoms  relaxed  towards  their  equilibrium 
distribution,  as  shown  in  fig.  2. 

The  profiles  in  Fig.  3  were  obtained  from  vertical  cuts  through  absorption  images.  They  provide 
clear  evidence  of  anti  ferromagnetic  interaction.  The  spin  structure  is  consistent  with  the  corresponding 
spin  domain  diagram  in  fig.  1  a.  Overlapping  mr-±l  clouds  as  observed  are  incompatible  with  the 
assumption  of  ferromagnetic  interaction. 

The  strength  c  =  (50  ±  20)  Hz  of  the  anti- ferromagnetic  interaction  was  estimated  by  determining 
C6,  the  location  of  the  mp  =  0  to  the  mp  =  ±1  boundary,  and  by  plotting  p  =  /zE'c^  versus  the 
quadratic  Zeeman  shift  q  =  qB^  as  shown  in  Fig.  4.  With  n  =  (2.9  ±  0.5)  x  10^'^  cm*"^  the  difference 
between  the  scattering  lengths  can  be  determined  to  no  =  3Aa  =  (3.5  ±  1.5)nB  =  (0.19±  0.08)  nm 
where  ap  denotes  the  Bohr  radius.  This  result  is  in  rough  agreement  with  a  theoretical  calculation 
of  02  —  flo  =  (5.5  ±  0.5)flB''‘-  The  anti  ferromagnetic  interaction  energy  corresponds  to  2.5  nK  in 
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our  condensates.  Still,  the  magnetostatic  (ferromagnetic)  interaction  between  the  atomic  magnetic 
moments  is  about  ten  times  weaker.  It  is  interesting  to  note  that  the  optically  trapped  samples  in 
which  the  domains  were  observed  were  at  a  temperature  of  the  order  of  100  iiK,  far  larger  than  the 
anti  ferromagnetic  energy.  The  formation  of  spin  domains  occurs  only  in  a  Bose  Einstein  condensate. 

Fig.  3  c  shows  a  profile  of  the  density  distribution  for  a  cloud  at  Bq  =  20  mG  and  almost  canceled 
gradient  (B'  <  2mG/cm).  No  mp  =  0  region  can  be  identified.  The  cloud  was  prepared  with  a  small 
total  angular  momentum.  Due  to  the  almost  -zero  gradient  and  the  non  zero  angular  momentum  the 
cloud  corresponds  to  a  point  in  the  shaded  region  in  Fig.  1  a,  rather  than  a  vertical  line  with  no  offset  as 
discussed  before  with  finite  gradients  and  zero  angular  momentum.  The  different  widths  of  the  profiles 
are  probably  caused  by  residual  field  inhomogenities.  Fig.  3  c  demonstrates  the  complete  miscibility  of 
the  mp  =  ±1  components. 

For  a  homogenous  two  component  system  the  criterion  for  miscibility  (immiscibility)  is 
when  the  mean  field  energy  is  parametrized  as  (27rrr/.V)(n^aa  +  + 

2na»6fla6)*  Here,  and  aa,b  are  densities  and  scattering  lengths  for  the  components  a  and  6,  and  the 
scattering  length  Oab  characterizes  the  interactions  between  particles  a  and  6.  In  our  spinor  condensate 
with  mixtures  of  the  mp  ±  1  components,  we  have  a_i  =  a^i  =  a  +  Sa  and  U-i+i  =  a  —  An.  Thus 
An  >  0,  like  experimentally  observed,  implies  miscibility.  For  a  mixture  of  the  mp  =  1  and  mp  =  0 
components,  we  find  no  =  n.  n+x  =  n  +  An  and  no+i  =  n  +  An,  corresponding  to  immiscibility.  For 
the  “^‘Rb  experiments"'*’  it  is  not  clear  whether  the  two  components  are  miscible  or  overlap  only  in  a 
surface  region  due  to  kinetic  energy-^. 

In  conclusion,  Bose  Einstein  condensates  of  sodium  occupying  all  three  hyperfine  states  of  the  F  =  1 
ground  state  multiplet  were  optically  trapped  in  low  magnetic  fields.  The  hyperfine  states  are  coupled 
by  spin  exchange  processes,  resulting  in  the  formation  of  spin  domains.  We  developed  spin  domain 
diagrams  for  both  the  anti  ferromagnetic  and  the  ferromagnetic  case,  and  showed  that  sodium  has  anti 
ferromagnetic  interactions,  whereas  the  opposite  case  is  predicted  for  the  *^'Rb  F  =  1  spin  multiplet'^ 
All  regions  in  the  spin  domain  diagrams  are  accessible  with  our  experimental  technique  and  thus  any 
combination  of  the  three  hyperfine  components  can  be  realized  by  applying  small  external  magnetic 
fields.  Of  special  interest  for  future  work  is  the  zero  magnetic  field  case,  where  the  rotational  symmetry 
should  be  spontaneously  broken.  We  observed  both  mi.scibility  and  immiscibility  of  hyperfine  compo¬ 
nents.  Thus  the  dynamics  and  possible  metastable  configurations*  of  two  interpenetrating,  miscible 
superfluid  components  with  arbitrary  admixtures  of  an  immiscible  component  (mp  =  0)  can 

now  be  studied. 
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Figure  captions: 


Figure  1:  Spin  domain  diagrams  for  spin  one  condensates.  The  structure  of  the  ground  state  spinor 
is  shown  as  a  function  of  the  linear  (-«  />)  and  quadratic  q)  Zeeman  energies.  Hyperfine  components 
are  mixed  inside  the  shaded  regions.  Solid  lines  indicate  a  discontinuous  change  of  state  populations 
whereas  dashed  lines  indicate  a  gradual  change.  The  behaviour  for  g  <  0  is  also  shown  although  it 
is  not  relevant  for  this  experiment.  For  c  =  0,  the  Zeeman  energy  causes  the  cloud  to  separate  into 
three  domains  with  mp  =  +1,0, —1  and  with  boundaries  at  \p\  =  g,  as  shown  in  b.  For  Co  ^  0,  the 
mean  field  energy  shifts  the  boundary  region  between  domains  and  leads  to  regions  of  overlapping  spin 
components.  In  the  anti  ferromagnetic  case  (a),  the  m/:'  =  0  component  and  the  wp=±\  components 
are  immiscible  (including  the  kinetic  energy  terms  in  Eqn.  (1)  would  lead  to  a  thin  boundary  layer)  and 
the  boundary  occurs  at  \p\  =  ^  +  c.  For  small  bias  fields,  with  q  <  c  and  |/>|  <  2c,  the  m/r  =  0  domain  is 
bordered  by  domains  in  which  components  are  mixed.  The  ratio  of  the  nip  =  ±l  populations 

in  these  regions  does  not  depend  on  q,  but  is  given  by  |4+|'/|^_|"  =  (2c  +  />)/(2c  —  p).  In  this  region  of 
small  fields,  the  boundary  to  the  m/:'  =  0  component  lies  at  |p|  =  2^^.  In  the  ferromagnetic  case  (c)  all 
three  components  are  generally  miscible,  and  have  no  sharp  boundaries.  Pure  /7iir  =  0  domains  occur  for 
|;v|  <  \/q(q  4|c|)  and  pure  nip  =  ±\  domains  for  |/>|  >  q.  Here,  in  contrast  to  the  anti  ferromagnetic 

case,  a  pure  nip  —  Q  condensate  is  skirted  by  regions  where  it  is  mixed  predominantly  with  either  the 
or  the  =  component.  The  contribution  of  the  third  component  is  very  small  (<  2%). 
In  all  mixed  regions  the  nif'zrO  component  is  never  the  least  populated  of  the  three  spin  components. 
This  qualitative  feature  can  be  used  to  rule  out  that  F  =  1  sodium  atoms  have  ferromagnetic  interactions. 


Figure  2:  Formation  of  ground  state  spin  domains.  Absorption  images  of  ballistically  expanding 
spinor  condensates  show  both  the  spatial  and  hyperfine  distributions.  Arbitrary  populations  of  the  three 
hyperfine  states  were  prepared  using  rf  transitions  (Landau  Zener  sweeps)  At  a  bias  field  of  about 
40  G  the  transitions  from  nip  =  —1  to  nip  =  0  and  from  nip  =  0  to  nip  =  +1  differ  in  frequency  by 
about  0.9  MHz  due  to  the  quadratic  Zeeman  shift  and  they  could  be  driven  separately.  The  images  of 
clouds  with  various  dwell  times  in  the  trap  show  the  evolution  to  the  same  equilibrium  for  condensates 
prepared  in  either  a  puire  nip  —  0  state  (upper  row)  or  in  equally  populated  nip  =  ±1  states  (lower 
row).  Between  os  to  los  dwell  time,  the  distribution  did  not  significantly  change,  although  the  density 
decreased  due  to  three  body  recombination.  The  bias  field  during  the  dwell  time  was  Bq  =  20 mG  and 
the  field  gradient  was  B'  =  llmG/cm.  These  images  were  taken  after  the  optical  trap  was  suddenly 
switched  off  and  the  atoms  were  allowed  to  expand.  Due  to  the  large  aspect  ratio  (typncally  60),  the 
expansion  was  almost  purely  in  the  radial  directions.  All  the  mean  field  energy  was  released  after  less 
than  1  ms,  after  which  the  atoms  expanded  as  free  particles.  Thus,  a  magnetic  field  gradient,  which  was 
applied  after  5  ms  time  of  flight  to  yield  a  Stern  Gerlach  separation  of  the  cloud,  merely  translated  the 
three  spin  components  without  affecting  their  shapes.  In  this  manner,  the  single  time  of  flight  images 
provided  both  a  spatial  and  spin  state  description  of  the  trapped  cloud.  Indeed,  the  shapes  of  the  three 
clouds  fit  together  to  form  a  smooth  total  density  distribution.  After  a  total  time  of  flight  of  25  ms 
the  atoms  were  optically  pumped  into  the  F  =  2  hyperfine  state  and  observed  using  the  nip  = +2  to 
m^  =  +3  cycling  transition.  This  technique  assured  the  same  transition  strength  for  atoms  originating 
from  different  spin  states.  The  size  of  the  field  of  view  for  a  single  spinor  condensate  is  1.7  mm  x  2.7  mm. 
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Figure  3:  Miscible  and  immiscible  spin  domains.  Axial  column  density  profiles  of  spinor  Bose> 
Einstein  condensates  are  shown,  obtained  from  time-of-flight  absorption  images  as  in  Fig.  *2.  The 
profiles  of  the  mr  =  ±1  components  were  shifted  to  undo  the  Stern-Gerlach  separation.  At  low  bias 
fields  (Fig.  3  a),  the  mp  =  0  component  was  skirted  on  both  sides  by  mp  =  ±1  components  with 
significant  n}p  =  ^^l  admixtures  thus  demonstrating  the  anti  ferromagnetic  interaction  (also  visible  in 
Fig.  2).  At  higher  fields  (Fig.  3  b),  the  mp  =  ±l  components  are  pushed  apart  further  by  a  larger 
wp  0  component  and  the  wp  =  admixtures  vanish.  They  could  not  be  resolved  for  quadratic 
Zeeman  energies  g  >  20  Hz.  The  anti  ferromagnetic  interaction  leads  to  immiscibility  of  the  iiip=^0  and 
the  7i}p  =  ±l  components.  The  kinetic  energy  in  this  boundary  region,  which  is  small  compared  to  the 
total  mean  field  energv,  is  released  in  the  axial  direction.  Due  to  this  axial  expansion  of  the  cloud  in 
the  time  of  flight  and  due  to  imperfections  in  the  imaging  system  including  the  limited  pixel  resolution, 
the  jnp  =  0  to  inp  =  ±1  boundary  is  not  sharp.  Fig.  3  c  demonstrates  the  complete  miscibility  of 
the  tup  =  ±1  components.  The  magnetic  field  parameters  were  Bq  =  20mG,  B'  =  llniG/cm  in  a, 
Bo  =  100  niG,  B'  =  11  mG/cm  in  b,  and  Bo  =  20mG,  B'  <  2mG/cm  in  c. 


Figure  4:  Estimate  of  the  anti  ferromagnetic  interaction  energy  c.  Plotted  is  the  linear  Zeeman 
energy  //  =  \fiB'z(,\  at  the  boundary  between  the  nip  =0  and  tiip  =  ±1  regions  versus  the  quadratic 
Zeeman  shift  q  =  qBq.  q  =  =  278Hz/G‘'.  and  pi  =  OsPib!"^^^  —  700kHz/G.  Qs  denotes 

the  electron  g  factor,  Uhjs  the  hyperfine  splitting  frequency,  and  pip  the  Bohr  magneton.  The  solid  line 
is  a  fit  of  the  function  Ip]  =  2^  ioi  q  <  c  and  lp|  =  ^  +  c  for  ^  >  c.  Extrapolating  the  linear  part 
to  zero  bias  field  (dashed  line)  yields  c  =  (50  ±  20)  Hz.  The  data  points  at  a  given  bias  field  represent 
p  =  pB'ct  for  different  gradient  fields  B'  and  thus  m^^O  regions  of  different  size.  The  scatter  of  these 
points  is  mainly  due  to  a  residual  magnetic  field  inhomogenities.  resulting  in  small  deviations  of  the 
local  gradient  B^  The  error  bar  represents  the  relative  error  of  all  data  points  of  30  %  in  p  and  5  %  in 
q  as  estimated  from  the  uncertainties  in  the  magnetic  field  calibration.  Furthermore  the  limited  pixel 
resolution  and  contributions  of  the  kinetic  energy  in  the  condensate  to  the  axial  expansion  enhance  the 
errors  for  the  determination  of  Zb  of  small  inp  =  0  regions. 
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We  present  a  method  of  adiabatically  changing  the  local  phase-space  density  of  an  ultracold  gas, 
using  a  combination  of  magnetic  and  optical  forces.  Applying  this  method,  we  observe  phase-space 
density  increases  in  a  gas  of  sodium  atoms  by  as  much  as  50-fold.  The  transition  to  Bose-Einstein 
condensation  was  crossed  reversibly,  attaining  condensate  fractions  of  up  to  30%.  Measurements  of 
the  condensate  fraction  reveal  its  reduction  due  to  interactions. 


The  physical  properties  of  atomic  gases  change  dra¬ 
matically  when  quantum  degeneracy  is  reached,  i.e.  when 
the  ground  state  population  approaches  unity  [1].  Re¬ 
cent  successes  in  reaching  quantum  degeneracy  with 
Bose  gases  [2-4]  have  relied  on  non-adiabatic,  irreversible 
methods  such  as  laser  and  evaporative  cooling.  The  pos¬ 
sibility  of  changing  the  ground  state  population  by  an 
adiabatic  (i.e.  isentropic)  change  in  the  trapping  poten¬ 
tial  had  been  overlooked  for  quite  some  time  [5].  In¬ 
deed,  in  the  case  of  an  ideal  gas,  adiabatic  changes  in 
the  strength  of  the  trapping  potential  do  not  change  the 
ground  state  population  [6,7].  However,  Pinkse  and  col¬ 
laborators  [8]  recently  showed,  both  theoretically  and 
experimentally,  that  by  changing  the  form  of  the  trap¬ 
ping  potential,  the  population  in  the  ground  state  can 
be  changed  without  changing  the  entropy.  For  a  non¬ 
degenerate  gas  the  ground  state  population  is  identical 
to  the  phase-space  density  F  =  where  n  is  the  den¬ 
sity  of  the  gas  and  Xp  is  the  thermal  de  Broglie  wave¬ 
length.  Within  the  type  of  trap  deformations  considered 
in  Ref.  [8]  the  maximum  increase  of  phase-space  density 
is  limited  to  a  factor  of  20. 

In  this  Letter,  we  show  that  a  more  general  deforma¬ 
tion  of  the  trapping  potential  can  increase  the  phase- 
space  density  by  an  arbitrary  factor,  and  implement 
this  scheme  using  a  combination  of  magnetic  and  optical 
forces.  Furthermore,  we  demonstrate  the  ability  to  cross 
the  Bose-Einstein  condensation  (BEC)  phase  transition 
reversibly. 

Adiabatic  increase  in  phase-space  density.  The 
type  of  trap  deformations  which  we  study  can  be  un¬ 
derstood  with  the  following  “two-box”  model.  Consider 
a  classical  gas  of  N  atoms  confined  in  a  box  of  volume 
Vq  zz  Vi  +  Vo  with  an  initial  phase-space  density  Fq.  Sup¬ 
pose  that  the  potential  within  a  sub- volume  Vo  of  the  box 
is  lowered  to  a  final  well-depth  U.  In  this  final  potential, 
the  gas  equilibrates  at  a  temperature  Tj,  and  the  density 
in  Vo  will  be  higher  than  that  in  Vi  by  the  Boltzmann 
factor  .  Using  the  condition  of  constant  entropy 

and  constant  particle  number,  one  obtains  the  relative 
increase  of  phase-space  density  in  Vo  compared  to  that 
in  Vo  before  compression: 


ln(Fo/Fo) 


UlksTj 

l  +  (Vo/Vi)e^A-Br,  * 


(1) 


For  deep  potential  wells,  where  U/ksTj  ^  ln(Vi/Vo), 
there  is  no  increase  in  phase-space  density  since  all  of  the 
gas  becomes  confined  in  Vo,  and  the  adiabatic  deforma¬ 
tion  corresponds  simply  to  a  uniform  compression  of  the 
gas.  For  shallow  potential  wells  (U/kBTf  <C  ln(Vi/Vo)), 
the  phase-space  density  in  Vo  increases  as  As 

is  varied  between  these  limits,  the  phase-space  den¬ 
sity  increase  reaches  a  maximum  which  is  greater  than 
(Vi/Vo)^/".  Thus,  by  choosing  an  extreme  ratio  of  vol¬ 
umes  Vi/Vo,  an  arbitrarily  large  increase  in  phase-space 
density  Fo/Fq  is  possible. 

To  demonstrate  this  phase-space  density  increase  in 
a  gas  of  trapped  atoms,  a  narrow  potential  well  (anal¬ 
ogous  to  Vo)  was  added  to  a  broad  harmonic  potential 
(corresponding  to  Vi)  by  focusing  a  single  infrared  laser 
beam  at  the  center  of  a  magnetic  trap.  First  a  gas  of 
atomic  sodium  was  evaporatively  cooled  to  a  temper¬ 
ature  higher  than  the  BEC  phase  transition  tempera¬ 
ture  in  the  cloverleaf  magnetic  trap  [9].  The  number  of 
atoms  and  their  temperature  were  adjusted  by  varying 
the  final  radio  frequency  (rf)  used  in  the  rf-evaporative 
cooling  stage  [10].  Afterwards,  the  magnetically  trapped 
cloud  was  decompressed  by  slowly  reducing  the  currents 
in  the  magnetic  trapping  coils.  Time-of-flight  absorption 
imaging  was  used  to  characterize  the  cloud.  The  total 
number  of  atoms  N  was  determined  by  integrating  the 
column  density  across  the  cloud,  and  the  initial  temper¬ 
ature  Jo  was  determined  by  one-dimensional  Gaussian 
fits  to  the  wings  of  the  density  distribution.  From  these, 
we  determined  the  fugacity  r  of  the  gas  by  the  relation 
gs{z)  =  N{huj/kBTo)^,  and  then  its  phase-space  density 

oo 

by  To  =  g3/2(~),  where  gn{z)  =  Y,  [H]-  Here,  Q  is 

t=l 

the  geometric-mean  trapping  frequency  of  the  magnetic 
trap,  as  determined  by  in  situ  measurements  [12].  These 
phase-space  density  measurements  were  calibrated  with 
images  from  magnetically  trapped  clouds  at  the  phase 
transition. 

The  optical  setup  was  similar  to  that  used  in  Ref.  [13]. 
The  infrared  laser  power  was  gradually  ramped- up  from 
zero  to  a  power  Pc  at  which  the  onset  of  BEC  was  seen 
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in  time-of-flight  images  of  clouds  released  from  the  de~ 
formed  trap;  this  implied  Tj  =  ^3/2(1)  =  2.612  for 
the  final  phase-space  density.  These  time-of-flight  im¬ 
ages  allowed  us  to  determine  the  final  temperature  Tj. 
The  depth  of  the  optical  potential  well  was  given  by 
Uc/f^B  =  37/iK  Pc/'^^o  (/^ni-/mW),  where  ivq  is  the  l/e“ 
beam-waist  radius  at  the  focus.  The  ramp-up  time  was 
made  long  enough  to  ensure  that  the  trap  deformation 
was  adiabatic,  but  also  short  enough  to  minimize  heat¬ 
ing  and  trap  loss.  Ramp-up  times  of  up  to  10  s  were 
used. 


FIG,  1,  Phase-space  density  increase  to  reach  BEC  vs, 
normalized  well-depth  (a).  Various  trap  settings  were  used: 
u?  =  27r  X  100  Hz,  Wo  =  9  fim  (triangles);  w  =  2t7  x  100  Hz, 
Wo  =  18 (open  circles);  and  u)  =  2;:  x  33  Hz,  wo  =  18 
(closed  circles).  The  solid  line  gives  the  prediction  of  Eq. 
(3).  (b)  Prior  to  deformation,  the  harmonic  trapping  poten¬ 
tial  V(r)  holds  a  cloud  above  the  BEC  transition  temperature 
(p  <  0).  (c)  When  a  potential  well  with  depth  Uc  =  is 
added,  a  small  condensate  forms. 


Fig.  la  shows  the  increases  in  phase-space  density 
which  were  measured  at  three  different  settings  of  the 
trap  parameters.  A  maximum  increase  by  a  factor  of  50 
was  obtained.  Condensates  were  observed  in  clouds  with 
temperatures  as  high  as  5  ^K.  Further  increases  were  hin¬ 
dered  by  limitations  in  laser  power  and  by  limits  to  the 
ramp-up  time  set  by  the  various  heating  and  loss  pro¬ 
cesses  in  the  deformed  trap.  The  scatter  in  the  data  is 
primarily  due  to  statistical  errors  at  the  level  of  30%  in 
our  measurements  of  Uc/^^bP/,  <^nd  in  the  determination 
of  the  transition  point. 

The  well-depth  Uc  required  to  reach  BEC  can  be  un¬ 
derstood  by  a  simple  model  depicted  in  Fig.  1.  We  begin 
with  a  harmonically  trapped  gas  with  a  chemical  poten¬ 
tial  fjt  <  0  (Fig.  lb).  Then,  by  lowering  a  potential  well 
to  a  depth  Uc  at  which  BEC  occurs,  one  essentially  mea¬ 
sures  n  diS  n  ~  —Uc  (Fig.  Ic).  However,  this  simple 
picture  neglects  the  change  in  the  chemical  potential  and 
the  temperature  of  the  gas  during  the  adiabatic  com¬ 
pression.  Thus,  let  us  consider  instead  the  condition  of 
constant  entropy.  The  entropy  per  non-condensed  par¬ 
ticle  in  a  harmonically  confined  Bose  gas  is  determined 
uniquely  by  its  fugacity  r  [8]: 
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This  equation  describes  the  entropy  of  the  gas  before 


compression,  with  the  fugacity  given  by  cq.  After  com¬ 
pression,  because  of  the  small  volume  of  the  poten¬ 
tial  well,  the  entropy  per  particle  is  approximately  that 
of  a  harmonically  trapped  gas  (Eq.  (2))  with  fugacity 
zj  =  Here  Tj  is  the  final  temperature  of 

the  gas.  Constant  entropy  then  implies  zo  =  zj.  The  ini¬ 
tial  phase-space  density  of  the  gas  is  Fq  =  ^3/2(‘^o)  and 
its  final  phase-space  density  is  F/  =  53/2(1)-  Thus  we 
obtain 


£/  =  g3/2(l) 

To  93/2{eM-Uc/kBTf))  • 

This  prediction,  shown  in  Fig.  la,  describes  our  data  well, 
and  accounts  for  the  universal  behaviour  of  our  measure¬ 
ments  over  a  wide  range  of  temperatures  and  well-depths. 

One  may  also  consider  the  process  of  adiabatically  in¬ 
creasing  the  phase-space  density  as  a  change  in  the  den¬ 
sity  of  states  D(e)  of  the  system.  By  increasing  the  well- 
depth  in  a  small  region  of  the  trap,  we  lowered  only  the 
energy  of  the  ground  state  and  a  few  excited  states.  Thus 
r,  a  local  quantity,  increases  as  the  ground  state  energy 
is  brought  closer  to  the  chemical  potential,  while  the  en¬ 
tropy,  a  global  property  of  the  gas,  is  unchanged  by  the 
minimal  modification  of  Z)(e). 


FIG.  2.  Momentum  distributions  of  the  thermal  cloud  at 
the  BEC  phase  transition  show  no  Bose-enhancement  for  the 
deformed  trap  (a),  but  a  clear  Bose-enhancement  for  the 
purely  magnetic  trap  (b).  Both  distributions  show  a  small 
condensate  peak.  Lines  show  a  Maxwell-Boltzmann  distri¬ 
bution  (dotted)  and  a  Bose- Einstein  distribution  (for  =  1) 
(dashed)  for  clouds  in  a  harmonic  trap,  which  were  fit  to  the 
thermal  wings.  The  momentum  distribution  is  shown  as  a 
profile  across  absorption  images  taken  after  40  ms  of  ballistic 
expansion. 

The  fact  that  global  properties  of  the  gas  are  not  af¬ 
fected  by  the  trap  deformation  can  also  be  seen  in  the 
momentum  distributions  probed  by  time-of-flight  imag¬ 
ing  [9,16].  The  onset  of  BEC  in  the  combined  opti¬ 
cal  and  magnetic  trap  is  signaled  only  by  the  formation 
of  a  condensate  peak.  The  remaining  thermal  cloud  is 
well  fit  by  a  Maxwell-Boltzmann  distribution,  which  de¬ 
scribes  a  magnetically  trapped  cloud  far  from  condensa¬ 
tion  (Fig.  2a).  In  contrast,  at  the  BEC  transition  in  the 
harmonic  magnetic  trap,  the  momentum  distribution  of 
the  thermal  cloud  is  clearly  Bose-enhanced  at  low  mo¬ 
menta  (Fig.  2b). 
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Adiabatic  condensation.  We  now  turn  to  the  stud¬ 
ies  of  adiabatic,  i.e.,  reversible,  condensate  formation.  A 
cloud  of  about  50  x  10^  atoms  was  cooled  to  the  the 
transition  temperature  in  the  magnetic  trap,  at  trap  fre¬ 
quencies  of  u!r  =  27r  X  20  Hz  and  uJz  =  27r  x  13  Hz  in 
the  radial  and  axial  direction,  respectively.  The  power 
of  the  infrared  laser  beam  (of  radius  tco  =  20 //m),  was 
ramped  up  over  1  s  and  held  at  a  constant  power  for  a 
dwell  time  of  1.5  s.  Condensate  fractions  as  small  as  1% 
could  be  distinguished  from  the  normal  fraction  by  their 
anisotropic  expansion  in  time-of-flight  images  [2,3].  The 
condensate  number  A^o  was  determined  by  subtracting 
out  the  thermal  cloud  background  using  Gaussian  fits  to 
the  thermal  cloud  in  regions  where  the  condensate  was 
clearly  absent. 


FIG.  3.  Condensate  formation  by  adiabatic  trap  deforma¬ 
tion.  The  condensate  fraction  is  plotted  against  the  (nor¬ 
malized)  well-depth.  Lines  show  predictions  for  an  ideal  gas 
(solid)  and  for  an  interacting  gas  (dashed). 

As  shown  in  Fig.  3,  the  adiabatic  trap  deformation 
yielded  condensate  fractions  of  up  to  15%;  we  observed 
condensate  fractions  of  30%  with  different  trap  settings. 
By  varying  the  dwell  time,  we  confirmed  that  clouds 
for  which  U/ksTj  <  1.5  suffered  no  significant  losses  of 
condensate  number  due  to  three- body  decay  or  heating, 
whereas  those  points  with  higher  values  of  U/ksTj  were 
affected  by  such  losses. 

For  an  ideal  Bose  gas,  the  result  of  this  adiabatic 
change  can  be  understood  as  follows.  Before  compres¬ 
sion,  the  cloud  of  lY  particles  at  the  BEC  transition  has 
an  entropy  Si  given  by  Eq.  (2)  as  Si  =  N  x  4g^{l)/g3{l). 
After  compression,  the  situation  is  similar  to  that  indi¬ 
cated  in  Fig.  Ic,  i.e.  because  of  the  small  volume  of  the 
attractive  well,  the  cloud  is  well-described  as  a  harmoni¬ 
cally  trapped  gas  with  //  =  — £/.  Thus,  accounting  for  the 
fact  that  condensate  particles  carry  no  entropy  and  again 
using  Eq.  (2),  which  gives  the  entropy  per  non-condensed 
particle,  we  equate  the  entropy  before  and  after  compres¬ 
sion  and  obtain 

_ %4(l)/g3(l) _ 

However,  this  simple  prediction  does  not  describe  our 
findings  well.  The  theory  described  above  has  two  short¬ 
comings.  First,  the  approximation  of  using  Eq.  (2)  for 


the  deformed  trap  is  not  strictly  valid.  However,  cal¬ 
culations  which  accounted  for  the  true  shape  of  the  de¬ 
formed  potential  changed  the  prediction  of  Eq.  (4)  only 
for  U jksTj  >  1,  and  only  slightly  improved  the  fit  to  our 
data. 

A  second  shortcoming  is  the  neglect  of  interactions.  It 
has  been  shown  that  the  condensate  fraction  of  harmon¬ 
ically  trapped  Bose  gases  with  repulsive  interactions  is 
reduced  in  comparison  to  that  of  an  ideal  gas  [17,18].  To 
estimate  this  effect  in  our  system,  we  use  the  ‘‘semi-ideal” 
model  of  Ref.  [18].  The  thermal  cloud  is  described  as  an 
ideal  gas  for  which  the  chemical  potential  is  raised  by 
griQ  =  47r7i^a7?o/m,  where  no  is  the  maximum  conden¬ 
sate  density,  m  the  mass  of  sodium,  and  a  =  2.75  nm  its 
scattering  length  [19].  This  simply  corresponds  to  using 
Eq.  (4)  with  the  substitution  U  ^  U  —  gno.  We  de¬ 
termined  no  using  the  Gross-Pitaevskii  equation  in  the 
Thomas- Fermi  limit  [20],  and  a  harmonic  approximation 
for  the  deformed  trap  potential  at  its  center.  Both  ap¬ 
proximations  are  valid  for  all  data  in  Fig.  3. 

This  approach  predicts  a  significant  reduction  of  the 
condensate  fraction  (Fig.  3,  dashed  line),  and  the  im¬ 
proved  agreement  with  our  data  is  strong  evidence  for 
this  effect.  In  contrast  to  related  studies  in  purely  har¬ 
monic  traps  [9,16],  which  did  not  show  evidence  for  inter¬ 
action  effects,  this  depletion  is  strongly  enhanced  by  the 
shape  of  the  deformed  potential.  The  mean-field  energy 
of  the  condensate  gno  is  large  because  the  condensate 
forms  in  the  tight  optical  potential,  while  the  transition 
temperature  Tc  is  small  since  it  is  determined  by  the  weak 
magnetic  potential.  The  remaining  discrepancy  between 
the  “semi-ideal”  prediction  and  our  data  may  be  due  to 
our  approximation  to  the  entropy  of  the  interacting  gas, 
which  fails  for  large  values  of  gno/ksTj  [21]. 

The  reversibility  of  crossing  the  BEC  phase  transition 
was  demonstrated  by  preparing  a  magnetically  trapped 
cloud  just  above  Tc.  We  then  sinusoidally  modulated 
the  power  of  the  infrared  light  at  1  Hz,  between  zero 
and  7  mW.  This  modulation  frequency  was  significantly 
smaller  than  the  magnetic  trap  frequencies  (u;^  =  27r  x 
48  Hz  and  u;.  =  27r  x  16  Hz).  These  low  frequencies  and  a 
large  optical  focus  (u’o  =  18  fim)  were  used  to  minimize 
trap  loss  due  to  inelastic  collisions. 

During  the  first  seven  condensation  cycles  the  conden¬ 
sate  fraction  oscillated  between  zero  and  6%i  (Fig.  4); 
later  probing  showed  repeated  condensation  for  at  least 
15  cycles.  The  peak  of  these  oscillations  decreased  slowly 
in  time.  The  temperature  oscillated  with  an  amplitude  of 
about  100  nK,  while  gradually  rising  by  about  10  nK/s, 
This  heating  and  the  decrease  of  the  peak  condensate 
fraction  are  consistent  with  similar  behaviour  in  clouds 
held  at  a  constant  infrared  power  [13].  Thus,  within  the 
stability  limitations  of  our  optical  setup,  the  repeated 
crossing  of  the  BEC  phase  transition  appears  fully  adia¬ 
batic. 
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FIG.  4.  Adiabatic  cycling  through  the  phase  transition. 
Shown  is  the  condensate  fraction  (a)  and  the  temperature 
(b)  vs  time  for  the  case  of  a  modulated  infrared  beam  (closed 
circles),  an  infrared  beam  ramped  up  to  a  constant  power 
(open  triangles),  and  no  infrared  light  (open  squares).  The 
solid  lines  are  guides  to  the  eye. 

This  method  of  creating  condensates  provides  insight 
into  their  formation,  which  was  recently  studied  experi¬ 
mentally  [22]  and  theoretically  [23].  For  example,  in  the 
experiment  described  above,  the  condensate  fraction  was 
found  to  lag  about  70  ms  behind  the  modulation  of  the 
laser  power,  which  is  a  measure  for  the  formation  time. 
That  this  lag  time  is  much  shorter  than  the  timescale 
for  deforming  the  potential  1  s)  confirms  that  the 
repeated  condensation  is  performed  adiabatically. 

In  other  experiments,  by  switching  on  the  infrared  light 
instantly,  we  observed  condensation  on  time  scales  much 
faster  than  the  oscillation  periods  in  the  magnetic  trap 
and  along  the  weakly  confining  axis  of  the  optical  trap. 
The  resulting  condensates  showed  striations  in  time-of- 
flight  images,  indicating  that  the  condensates  formed  into 
excited  states  of  the  deformed  potential.  Such  studies  of 
shock-condensation  might  give  new  insight  into  the  for¬ 
mation  of  quasi-condensates  and  condensation  into  ex¬ 
cited  states  [24,25]. 

In  conclusion  we  have  demonstrated  the  adiabatic 
Bo.se- Einstein  condensation  of  an  ultracold  gas  of  atomic 
sodium.  Changes  in  the  trapping  potential  resulted  in 
large  phase-space  density  increases  and  allowed  for  re¬ 
peated  crossings  of  the  BEC  phase  transition.  This 
method  allows  for  detailed  studies  of  condensate  forma¬ 
tion  and  the  phase  transition.  The  combined  trapping 
potential  widens  the  range  of  trap  parameters  over  which 
BEC  can  be  studied.  This  was  used  to  strikingly  enhance 
the  role  of  interactions,  and  led  to  higher  transition  tem¬ 
peratures  (up  to  5  fiK)  than  achieved  in  purely  magnetic 
traps. 
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Present  experiments  with  Bose  condensed  gases  can  be  largely  described  by  a  semi-ideal  two-gas 
model.  In  this  model,  the  condensate  is  influenced  only  by  the  mean-field  repulsion  among  condensed 
atoms,  while  the  thermal  cloud  is  considered  an  ideal  gas  confined  by  cin  effective  potential  that 
consists  of  the  external  trap  and  the  mean- field  repulsion  by  the  condensate.  This  simple,  intuitive 
model  provides  explicit  analytical  expressions  for  the  density  distributions  of  the  condensate  and 
the  thermal  component.  It  describes  the  reduction  of  the  condensate  fraction  relative  to  that  of  an 
ideal  gas  as  a  consequence  of  the  positive  chemical  potential  due  to  interactions  in  the  condensate. 
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Recent  experiments  with  Bose-Einstein  condensed  al¬ 
kali  gases  (cf.  e.g.  [1,  2])  have  been  successfully  de¬ 
scribed  by  different  approximation  schemes.  The  dilute¬ 
ness  of  these  gases,  expressed  as  na^  1  where  a  is  the 
two-body  scattering  length,  allows  a  description  in  terms 
of  an  effective  mean-field  theory.  At  zero  temperature 
the  condensate  can  be  described  by  the  Gross-Pitaevskii 
equation  [3,  4],  a  nonlinear  Schrodinger  equation.  In 
present  experiments  interactions  are  strong  enough  so 
that  the  chemical  potential  p  is  much  larger  than  the 
level  spacings  hujx,  Tiujyj  cind  fiujz  of  the  external  har¬ 
monic  potential.  Therefore  the  use  of  the  Thomas-Fermi 
approximation  [5,  6],  which  neglects  the  kinetic  energy  of 
the  condensate,  has  lead  to  a  tremendous  simplification 
in  the  understanding  of  these  experiments. 

Recently,  much  theoretical  and  experimental  atten¬ 
tion  has  been  focused  on  finite  temperature  properties 
of  Bose-condensed  gases.  A  semiclassical  Hartree-Fock- 
Popov  description  has  produced  excellent  agreement  with 
experimentally  measured  equilibrium  quantities  [7].  Fur¬ 
ther,  unlike  in  a  spatially  homogeneous  gas,  it  was  found 
that  collective  elementary  excitations  contribute  very  lit¬ 
tle  to  thermodynamic  properties  [7]  due  to  the  relatively 
small  spatial  overlap  between  condensate  and  thermal 
component  in  a  harmonic  trap.  Thus,  even  a  semiclas¬ 
sical  Hartree-Fock  approximation  is  a  fairly  accurate  de¬ 
scription  of  these  gases,  except  for  a  small  temperature 
range  close  to  the  critical  temperature.  For  example,  ex¬ 
cellent  agreement  has  been  achieved  by  a  comparison  of 
Hartree-Fock  density  profiles  with  the  results  of  a  com¬ 
prehensive  Path  Integral  Monte  Carlo  calculation  [8]. 

In  this  paper  we  introduce  a  further  simplification 
of  the  self-consistent  Hartree-Fock  model  by  neglecting 
atomic  interactions  in  the  thermal  component  of  the  gas. 
Like  the  introduction  of  the  Thomas-Fermi  approxima¬ 
tion  [5]  the  reduction  to  Eqs.  (4)-(6)  and  (12)  has  the 
important  advantage  of  providing  an  explicit  analytical 
description  of  a  trapped  Bose  gas  at  finite  temperature. 


We  consider  atoms  that  are  trapped  in  a  generally 
anisotropic  harmonic  potential.  Under  the  appropriate 
rescaling  of  energy  and  lengths,  the  harmonic  potential 
has  the  form  V(r)  =  r-/2.  Here,  the  dimensionless  carte¬ 
sian  coordinates  of  r  are  given  in  units  of  yjhuj/imvl  y  .,, 

with  uj  =  Other  lengths,  like  the  thermal 

wavelength  Xp  and  the  scattering  length  a,  are  scaled  in 
units  of  the  natural  length  \Jh/inuj  of  a  harmonic  oscil¬ 
lator  with  angular  frequency  uj.  The  energy  V{r)  as  well 
as  other  energies  are  implicitly  given  in  units  of  hij. 

Then,  in  the  self-consistent  Hartree-Fock  model,  the 
densities  of  the  condensate  and  the  thermal  component 
are  given  as  solutions  of  the  two  coupled  equations  [5,  7] 

„o(r)  =  M  -  ry2  -  •2C^nT(r)  _  2UnT(r))  (1) 

nrCr)  =  ^53/2  .  (2) 

Here,  the  Thomas-Fermi  approximation  has  been  applied 
which  neglects  the  kinetic  energy  of  the  condensate  and 
implicitly  finite  size  effects.  The  strength  of  the  atomic 
interactions  is  given  by  U  ^  Attg.  We  use  the  thermal 
wavelength  Xp  =  and  the  Bose  function  is 

defined  by  ga{^)  =  YlJLi  chemical  potential 

p  is  determined  by  the  constraint  of  a  given  total  atom 
number 

N  =  j dr  [no(r)  +  nT(r)].  (3) 

For  further  reference,  Eqs.  (1)  -  (3)  will  be  denoted  as 
the  interacting  model.  Even  though  it  represents  already 
a  relatively  concise  description,  these  equations  still  have 
to  be  solved  self-consistently.  This  is  done  numerically 
by  an  iterative  procedure  [7] . 

If  one  further  neglects  the  mean-field  repulsion  from 
non-condensed  atoms  Eqs.  (1)  and  (2)  are  solved  by  the 
explicit  relations 


1 


no(r)  =  ^^-^^(//-rV2)  (4) 

«T(r)=^53/2(e-l’-'/=^-'‘'/'‘"^),  (5) 

with  the  chemical  potential  (i  given  by 

,6) 

These  relations  summarize  a  simple  model  of  the  partly 
condensed  gas  which  we  denote  as  the  semi-ideal  model. 
The  condensate  density  no  is  that  of  a  zero  temperature 
condensate  with  Nq  atoms,  while  the  density  ut  of  the 
non-condensed  cloud  is  that  of  an  ideal  gas  of  bosons 
confined  in  the  combination  of  the  external  potential  and 
the  repulsive  mean-field  potential  due  to  the  condensate 
atoms.  Its  wide  range  of  validity  has  been  confirmed  by 
a  numerical  analysis  in  Ref.  [9].  However,  in  order  to 
make  it  an  intuitively  appealing,  usable  description  of 
a  trapped  Bose  gas,  the  condensate  fraction  has  to  be 
determined  analytically  as  a  function  of  temperature. 

For  later  convenience  let  us  introduce  the  reduced 
chemical  potential  fi 


Here,  the  critical  temperature  Tc  is  that  of  a  trapped 
ideal  gas  [10,  11]  keTc  =  where  ^^(l)  =  C(q) 

was  expressed  in  terms  of  the  Riemann  Zeta  function. 
The  reduced  chemical  potential  depends  on  the  scaling 
parameter  t]  [7] 

»?  =  gfr  =  ic(3)'/"(l5iV'^®a)'^'  (8) 

which  describes  the  strength  of  the  atomic  interactions 
within  the  condensate.  The  scaling  parameter  is  indepen¬ 
dent  of  the  system  size  when  the  thermodynamic  limit  is 
taken  in  the  usual  way  [N  oo,a;  — )■  0,  Nuj^  =  const.). 
Due  to  its  relatively  weak  dependence,  rj  assumes  a 
value  close  to  0.3  in  most  recent  experiments  [7]. 

The  condensate  fraction  is  determined  by  integration 
over  the  thermal  density  distribution.  For  T  <  the  re¬ 
sult  can  be  written  in  terms  of  incomplete  Gamma  func¬ 
tions  [12] 


N  '  v^C(3) 


LJ-fLy 
^c(3)  \tJ 

noo 

+  /  dev/e 

J  ^ 


{c-n)lkBT 


) 

+  e^'^r(3/2,j//)]. 


(9) 


Eq.  (9)  is  still  an  implicit  expression  since  it  depends 
again  on  the  reduced  chemical  potential.  Therefore,  fi 
and  Nq/N  are  given  as  simultaneous  solutions  of  Eqs.  (7) 
and  (9).  For  temperatures  T  >  Tc  condensate  does  not 
exist  and  the  reduced  chemical  potential  //  <  0  is  trivially 
defined  by  the  single  equation  g3{e^)  =  C(3)(T/Tc)“^. 

Let  us  consider  the  results  of  the  semi-ideal  model  in 
the  limits  of  small  and  large  //,  corresponding  to  T 
and  T  ^  0,  Eq.  (9)  is  therefore  rewritten 


A^o  ^  .  2  1 

N  v^C(3) 


(10) 


in  terms  of  the  function 

pX  poo 

f(x)  =  j  +  J  (11) 

where  the  reduced  energy  e  =  c/ksT  has  been  intro¬ 
duced.  The  function  can  be  expressed  as  a  power 

series  in  where  p  runs  from  0  (T  =  Tc)  to  oo 

(T  =  0).  Truncating  the  series  after  the  lowest  nontrivial 
order  \delds 


The  first  two  terms  of  Eq.  (12)  correspond  to  the  conden¬ 
sate  fraction  of  an  ideal  gas,  whereas  the  third  describes 
the  influence  of  the  condensate  repulsion  to  lowest  order 
in  p.  This  expansion  adequately  describes  the  solution  of 
Eq.  (9)  over  the  entire  range  of  temperatures,  provided 
that  T]  1,  A  further  simplification  is  derived  by  solving 
Eq.  (12)  to  lowest  order  in  arriving  at 


(13) 

Within  Eqs,  (12)  and  (13),  the  only  effect  of  the  con¬ 
densate  interactions  is  to  shift  the  chemical  potential  to 
positive  values,  raising  the  bottom  of  the  effective  poten¬ 
tial  which  the  thermal  atoms  experience.  The  effect  of 
changing  the  shape  of  the  potential  enters  only  at  higher 
orders  of  the  expansion.  Eq.  (12)  is  thus  a  good  ap¬ 
proximation  of  Eq.  (9)  if  most  of  the  thermal  atoms  are 
spatially  separated  from  the  condensate,  a  condition  sat¬ 
isfied  for  temperatures  close  to  Tc  and  for  small  rj, 

Eq.  (12)  is  not  expected  to  be  a  good  description  for 
extremely  low  temperatures.  The  understanding  of  this 
limit  is  facilitated  by  f{x)  obeying  the  differential  equa¬ 
tion 


dx^ 


f  =  f-2s/^. 


(14) 
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The  function  f{x)  has  to  grow  slower  than  for  large 
values  of  x.  Otherwise,  the  power  series  in  Eq.  (10)  would 
not  converge.  As  a  consequence,  the  second  derivative  in 
Eq.  (14)  vanishes  for  large  x.  Thus,  the  asymptotic  limit 
of  the  condensate  fraction  for  low  temperatures  is  given 
by  f(x)  =  2y/x  or 


4  C(5/2) 
y/^  C{3) 


(15) 


However,  it  has  been  shown  [7]  that  the  leading  term  of 
the  uncondensed  fraction  around  zero  temperature  scales 
like 


(tV 

N  ''2v/2C(3)  \Tc) 


(16) 


due  to  quasi-particle  contributions  which  are  included  in 
neither  the  semi-ideal  nor  the  interacting  Hartree-Fock 
treatment. 

According  to  Fig.  1,  where  t]  =  0.31  has  been  assumed, 
the  high  temperature  expansions  Eqs.  (12)  and  to  some 
extent  Eq.  (13)  are  very  good  approximations  of  the  exact 
condensate  fraction  Eq.  (9)  of  the  semi-ideal  model  over 
the  whole  range  of  temperatures  T  <  Tc.  In  contrast, 
we  have  found  that  the  validity  of  the  low  temperature 
limits  of  Eqs.  (15)  and  (16)  are  restricted  to  a  very  small 
range  around  zero  temperature  which  makes  them  nearly 
irrelevant  for  practical  purposes. 

The  condensate  fraction  given  by  the  approximate  so¬ 
lution  of  the  semi-ideal  model  in  Eq.  (12)  is  compared 
in  Fig.  2  with  the  result  of  an  iterative  solution  of  the 
interacting  model  (Eqs.  (l)-(3)).  The  difference  between 
the  curves  is  small  except  near  the  critical  temperature, 
where  the  interactions  among  thermal  atoms  shift  the 
critical  temperature  [13].  This  shift  is  not  exhibited  by 
the  semi-ideal  model.  The  wide  range  of  validity  of  the 
semi-ideal  model  has  also  been  confirmed  by  a  numer¬ 
ical  comparison  with  a  Hartree-Fock-Popov  calculation 
[9].  The  latter  includes  collective  excitations,  which  are 
ignored  in  the  models  shown  in  Fig.  2.  These  lead  to 
an  additional  but  negligible  decrease  of  the  condensate 
fraction  [7]. 

A  direct  comparison  between  the  density  distributions 
of  the  different  models  is  given  in  Fig.  3  for  T/Tc  =  0.5. 
Here,  the  Thomas-Fermi  solution  of  the  condensate  den¬ 
sity  of  Eq.  (1)  has  been  replaced  by  a  solution  of  the  full 
Gross-Pitaevskii  equation.  The  main  implication  of  the 
included  kinetic  energy  of  the  condensate  is  a  smoothed 
out  condensate  surface.  Within  the  Thomas-Fermi  ap¬ 
proximation  almost  no  difference  is  seen  between  the  two 
curves.  Even  though  the  unaccounted  existence  of  collec¬ 
tive  excitations  changes  the  densities  of  condensate  and 
thermal  component  close  to  the  center  of  the  trap  by  a 
significant  amount,  their  influence  on  the  total  density  is 
negligibly  small  for  most  temperatures  [8]. 

In  recent  experiments,  trapped  Bose  gases  have  been 
probed  in-situ  using  either  non-destructive  dispersive 


imaging  techniques  [14],  or  off-resonant  absorption  imag¬ 
ing  [15].  These  methods  are  important  new  tools  in  the 
understanding  of  trapped  Bose  gases,  with  distinct  ad¬ 
vantages  over  the  previously  used  time  of  flight  measure¬ 
ments  of  expanding  atom  clouds  [1,  2].  We  therefore 
expect  that  the  results  of  this  paper  will  become  useful 
for  future  in-situ  experiments  with  trapped  Bose  gases. 
These  new  optical  detection  techniques  measure  column 
or  line  density  profiles;  we  therefore  determine  these  pro¬ 
files  from  Eqs.  (4)  and  (5)  by  integration  over  one  or  two 
dimensions. 

In  cylindrical  coordinates,  where  the  probe  beam  is 
directed  along  the  z-axis,  the  according  column  densities 
are  given  by 


no{p) 


4v/2  {p  - 


nrip)  = 


u 


0{p-p-/-2) 


C  <ii;^93/2(e  I' 


p{p)  <  0 
p{p)  >  0 


(17) 

(18) 


where  fi{p)  =  (//  —  p^/2)/kBT.  Equivalently,  we  obtain 
the  line  densities 


no(z)  =  J^  £/?]:  _  ..2/2) 

„  -  (^^)~  /  55/2(e'‘‘<"^)  pi:)  <  0 

A5  \2C(5/2)-55/2(e-^(^))  /l(--)>0 


(19) 

(20) 


with  p{z)  =  (p-  zy2)/kBT. 

We  therefore  conclude  that  the  proposed  semi-ideal 
model  of  Eqs.  (4)“(6)  and  (12)  represents  a  rather  good 
description  of  a  trapped  Bose  gas  at  finite  temperature 
provided  that  p  is  not  larger  than  the  values  rj  ^  0.3 
to  0.4  achieved  in  recent  experiments.  Its  mathematical 
simplicity  gives  a  much  clearer  conceptual  picture  than 
the  self-consistent  approaches  that  have  been  used  be¬ 
fore.  The  observed  agreement  of  the  semi-ideal  model 
with  a  self-consistent  Hartree-Fock  calculation  has  con¬ 
firmed  that  the  density  of  the  thermal  component  is  too 
low  to  necessitate  an  interacting  gas  description.  Instead, 
the  observed  depletion  of  the  condensate  fraction  arises 
from  condensate  interactions,  leading  to  a  positive  value 
of  the  chemical  potential  or  equivalently  to  a  raised  bot¬ 
tom  of  the  effective  potential  in  which  the  thermal  com¬ 
ponent  resides.  It  is  well  described  by  Eq.  (12),  which  is  a 
low  order  expansion  in  ft.  In  contrast,  the  changed  shape 
of  the  effective  potential  does  not  influence  the  thermal 
population  in  a  noticeable  way. 
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FIG.  1.  Condensate  fraction  as  a  function  of  temperature, 
given  by  the  semi-idecil  model.  The  dashed  line  shows  the 
exact  solution  Eq.  (9)  of  the  semi-ideal  model.  The  solid  line 
represents  its  high-temperature  limit  Eq.  (12).  The  explicit 
expression  of  Eq.  (13)  is  given  by  the  dotted  line.  This  plot 
demonstrates  the  validity  of  Eq.  (12)  over  the  whole  range  of 
temperatures  T  <Tc.  The  explicit  approximation  of  Eq.  ( 13) 
fails  only  in  the  vicinity  of  the  critical  temperature  (shown  in 
the  inset). 
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FIG.  2.  Condensate  fraction  plotted  as  a  function  of  tem¬ 
perature,  comparing  the  discussed  models.  The  ideal  gas 
condensate  fraction  is  given  by  the  dashed-dotted  line.  The 
dotted  line  is  a  numerical  solution  of  the  Hartree-Fock  Eqs. 
(l)-(3).  The  solid  line  is  obtained  from  the  analytical  ap- 
pro.ximation  Eq.  (12),  corresponding  to  the  solid  line  of  Fig. 
1.  Except  near  the  critical  temperature  (shown  in  the  inset), 
the  semi-ideal  model  well  approximates  the  solution  of  the 
interacting  model.  The  critical  temperature  Tc  refers  to  the 
case  of  an  ideal  gas. 


FIG.  3.  A  comparison  of  the  density  distributions  at 
T/Tc  =  0.5  predicted  by  the  two  discussed  models.  The  up¬ 
per  solid  line  shows  the  density  of  the  thermal  component 
as  given  by  Eq,  (5).  In  comparison,  the  lower  solid  line  has 
been  calculated  numerically  using  Eq.  (2).  However,  a  nu¬ 
merical  solution  of  the  full  Gross-Pitaevskii  equation  has  been 
used  instead  of  the  Thomas- Fermi  approximation  of  Eq.  (1). 
A  small  difference  exists  at  the  edge  of  the  condensate  due 
to  the  neglect  of  the  kinetic  energy  in  the  semi-ideal  model. 
However,  the  total  density  is  only  weakly  affected  by  this  ef¬ 
fect,  as  is  seen  by  the  two  dcished-dotted  lines.  Again,  the 
upper  dashed-dotted  line  refers  to  the  total  density  given  by 
the  semi-ideal  model  while  the  lower  curve  corresponds  to 
the  self-consistent  Hartree-Fock  approximation.  Here,  a  total 
atom  number  of  =  5  x  10^  was  assumed. 
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I  Observation  of  Feshbach 

resonances  in  a 

Bose- Einstein  condensate 

S.  lnouye»,  M.  R.  Andrews^,  J.  Stenger,  H.-J.  Miesner»,  D.  M.  Stamper-Kurn*  &  W.  Ketterle* 

*  Department  of  Physics  and  Research  Laboratory  of  Bectronics,  Massachusetts  Institute  of  Technology,  Cambridge,  Massachusetts  02139,  USA 

It  has  long  been  predicted  that  the  scattering  of  ultracold  atoms  can  be  altered  significantly  through  a  so^alled 
‘Feshbach  resonance’.  "nMO  such  resonances  have  now  been  observed  in  optically  trapped  Bose-Einstein  condensates 
of  sodium  atoms  by  varying  an  external  magnetic  field.  They  gave  rise  to  enhanced  inelastic  processes  and  a  dispersive 
variation  of  the  scattering  length  by  a  factor  of  over  ten.  These  resonances  open  new  possibilities  for  the  study  and 
manipulation  of  Bose-Einstein  condensates. 


I 
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Bose-Einstein  condensates  of  atomic  gases  offer  new  opportunities 
for  studying  quantum-degenerate  fluids‘''\  All  the  essential  proper¬ 
ties  of  Bose  condensed  systems — the  formation  and  shape  of  the 
condensate,  the  nature  of  its  collective  excitations  and  statistical 
fluctuations,  and  the  formation  and  dynamics  of  solitons  and 
vortices — are  determined  by  the.  strength  of  the  atomic  interac¬ 
tions.  In  contrast  to  the  situation  for  superfluid  helium,  these 
interactions  are  weak,  allowing  the  phenomena  to  be  theoretically 
described  from  ‘first  principles’.  Furthermore,  in  atomic  gases  the 
interactions  can  be  altered,  for  instance  by  employing  different 
species,  changing  the  atomic  density,  or,  as  in  the  present  work, 
merely  by  varying  a  magnetic  field. 

At  low  temperatures,  the  interaction  energy  in  a  cloud  of  atoms  is 
proportional  to  the  density  and  a  single  atomic  parameter,  the 
scattering  length  a  which  depends  on  the  quantum -mechanical 
phase  shift  in  an  elastic  collision.  It  has  been  predicted  that  the 
scattering  length  can  be  modified  by  applying  external  magnetic^’'®, 
optical*^’*'  or  radio- frequency (r.f.)  fields.  Those  modifications  are 
only  pronounced  in  a  so-called  “Feshbach  resonance”’**,  when  a 
quasibound  molecular  state  has  nearly  zero  energy  and  couples 
resonantly  to  the  free  state  of  the  colliding  atoms.  In  a  time- 
dependent  picture,  the  two  atoms  are  transferred  to  the  quasibound 
state,  ‘stick’  together  and  then  return  to  an  unbound  state.  Such  a 
resonance  strongly  affects  the  scattering  length  (elastic  channel),  but 
also  affects  inelastic  processes  such  as  dipolar  relaxation*^  '  and  three- 
body  recombination.  Feshbach  resonances  have  so  far  been  studied 
at  much  higher  energies*^  by  varying  the  collision  energy,  but  here 
we  show  that  they  can  be  ‘tuned’  to  zero  energy  to  be  resonant  for 
ultracold  atoms.  The  different  magnetic  moments  of  the  free  and 
quasibound  states  allowed  us  to  tune  these  resonances  with  mag¬ 
netic  fields,  and  as  a  result,  minute  changes  in  the  magnetic  field 
strongly  affected  the  properties  of  a  macroscopic  system. 

Above  and  below  a  Feshbach  resonance,  the  scattering  length  a 
covers  the  full  continuum  of  positive  and  negative  values.  This 
should  allow  the  realization  of  condensates  over  a  wide  range  of 
interaction  strengths.  By  setting  one  can  create  a  condensate 
with  essentially  non-interacting  atoms,  and  by  setting  <i  <  0  one  can 
make  the  system  unstable  and  observe  its  collapse.  Rapid  tuning  of 
an  external  magnetic  field  around  a  Feshbach  resonance  will  lead  to 
sudden  changes  of  the  scattering  length.  This  opens  the  way  to 
studies  of  new  dynamical  effects  such  as  novel  forms  of  collective 
oscillations  or  the  sudden  collapse  of  a  large  condensate  when  the 
scattering  length  is  switched  from  positive  to  negative’^ 
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Theoretical  predictions 

Calculations  for  Feshbach  resonances  in  external  magnetic  fields 
have  been  reported  for  the  lower  hyperfine  states  of  the  atoms  Li 
(ref.  8),  K  (ref.  10),  Na  (ref.  8),  Rb  (ref.  9)  and  Cs  (refs  6, 7).  They  are 
typically  spaced  by  several  hundred  gauss,  and  for  Li  and  Na  occur 
outside  the  range  where  states  in  the  lower  hyperfine  manifold  are 
weak-field-seeking  and  can  be  magnetically  trapped.  Recent  experi¬ 
mental  efforts  to  observe  Feshbach  resonances  have  concentrated  on 
^”Rb  (ref.  17)  and  on  ®^Rb  (ref.  18  and  C.  E.  Wieman,  personal 
communication)  where  Feshbach  resonances  hav^e  been  predicted  at 
relatively  low  magnetic  fields^.  However,  our  recently  demonstrated 
all-optical  confinement  of  a  Bose  condensate*^  opened  the  possibi¬ 
lity  of  observing  Feshbach  resonances  for  strong- field -seeking  states  ' 
which  cannot  be  trapped  in  a  d.c.  magnetic  trap.  The  optical 
trapping  potential  is  unaffected  by  magnetic  fields  and  is  indepen-  | 
dent  of  the  hyperfine  ground  state.  We  report  here  the  observation  I 
of  two  Feshbach  resonances  of  sodium  in  a  strong-field-seeking 
state. 

Several  Feshbach  resonances  in  sodium  are  caused  by  quasibound  i 
h>"perfine  states  of  the  second  highest  vibrational  level,  v  =  14,  of 
the  triplet  potential  of  the  sodium  dimer.  The  lowest  magnetic  field 
value  Bq  for  a  strong  Feshbach  resonance  in  sodium  was  predicted  to 
lie  in  the  range  760  <  Bq  <  925  G  (B.  J.  Verhaar  and  F.  A.  van  j 
Abeelen,  personal  communication).  It  occurs  in  collisions  between  | 
atoms  in  the  lowest  hyperfine  state  |m5=  -  1/2.  m;  = +3/2>,  | 
which  correlates  with  the  (F  =  1,  tUf  =  -f  1)  state  at  low  fields  (S,  /  * 
and  f  are  the  usual  quantum  numbers  for  the  electronic,  nuclear  I 
and  total  spin,  respectively).  This  Feshbach  resonance  is  due  to  a  ; 
quasibound  molecular  state  |S=1,mi5  =  +!,/=!, =  +l).  A  | 
much  weaker  resonance  due  toa|5  =  =  tUf  =  -fl)  i 

state  (which  is  almost  degenerate  with  the  other  quasibound  state)  | 
was  predicted  to  occur  50  to  75  G  below.  | 

Near  a  Feshbach  resonance,  the  scattering  length  a  should  vary  I 
dispersively  as  a  function  of  magnetic  field  B  (ref.  8):  j 

-  =  “(‘■5^)  ‘"I 

where  id  parametrizes  the  width  of  the  resonance  at  B  =  B^,  and  a  is 
the  scattering  length  outside  the  resonance.  For  sodium,  a  was 
found  spectroscopically  to  be  2.75  nm  at  zero  field,  and  increases  to 
the  triplet  scattering  length  of  4.5  nm  (ref.  20)  at  high  magnetic 
fields.  The  widths  A  for  the  strong  and  weak  resonance  were  j 
predicted  to  be  1 G  and  0.01  G,  respectively  (B.  J.  Verhaar  and  F.  A.  | 
van  Abeelen,  personal  communication).  i 
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I  Experimental  set-up 

j  Bose-Einstein  condensates  in  the  |F  =  1,  Wp  =  —  l)  s^ate  were 
!  produced  as  in  our  previous  work  by  laser  cooling,  followed  by 
I  evaporative  cooling  in  a  magnetic  trap^^  The  condensates  were 
’  transferred  into  an  optical  dipole  trap  formed  at  the  focus  of  an 
infrared  laser  beam*^  Atoms  were  then  spin-flipped  with  nearly 
100%  efficiency  to  the  |f  =  1,  =  +l)  state  with  an  adiabatic  r.f 

sweep  while  applying  a  1  G  bias  field.  Without  large  modifications  of 
our  magnetic  trapping  coils,  we  could  provide  bias  fields  of  up  to 
j  — 1,200  G,  but  only  by  using  coils  producing  axial  curvature^^ 

I  which  for  high-field-seeking  states  generated  a  repulsive  axial 
!  potential.  At  the  highest  magnetic  fields,  this  repulsion  was  stronger 
j  than  the  axial  confinement  provided  by  the  optical  trap.  To  prevent 
.  I  the  atoms  from  escaping,  two  ‘end-caps’  of  far-off-resonant  blue- 
j  detuned  laser  light  were  placed  at  the  ends  of  the  condensate, 
^  creating  a  repulsive  potential,  and  confining  the  atoms  axially 
‘  j  (Fig.  la).  For  this,  green  light  at  514  nm  from  an  argon-ion  laser 
I  was  focused  into  two  sheets  about  200  jxm  apart.  The  focus  of  the 
;  optical  trap  was  placed  near  the  minimum  of  the  bias  field  in  order 
I  to  minimize  the  effect  of  the  destabilizing  magnetic  field  curvature. 

■  The  axial  trapping  potential  at  high  fields  was  approximately  “W”- 
:  shaped  (Fig.  lb),  and  had  a  minimum  near  one  of  the  end-caps  as 
I  obser\^ed  by  phase- contrast  imaging^  (Fig.  Ic,  d). 
j  The  calibration  factor  between  the  current  (up  to  -^400  A)  in  the 
;  coils  and  the  magnetic  bias  field  was  determined  with  an  accuracy  of 
!  2%  by  inducing  r.f.  transitions  within  the  \F  =  l)  ground-state 
h\'perfine  manifold  at  about  40  G.  Additionally,  an  optical  reso¬ 
nance  was  found  around  1,000  G,  where  the  Zeeman  shift  equalled 
;  the  probe  light  detuning  of  about  1.7  GHz  and  led  to  a  sign-reversal 
of  the  phase-contrast  signal.  These  two  calibrations  agreed  within 
their  uncertainties. 

The  condensate  was  observed  in  the  trap  directly  using  phase- 
contrast  imaging’"  or  by  using  time-of-flight  absorption 
imaging'  ""'.  In  the  latter  case,  the  optical  trap  was  suddenly 
switched  off,  and  the  magnetic  bias  field  was  shut  off  1-2  ms  later 
to  ensure  that  the  high-field  value  of  the  scattering  length  was 
responsible  for  the  acceleration  of  the  atoms.  After  ballistic  expan¬ 
sion  of  the  condensate  (either  12  or  20  ms),  the  atoms  were  optically 
pumped  into  the  |f  =  2)  ground  state  and  probed  using  resonant 
light  driving  the  cycling  transition.  The  disk-like  expansion  of  the 


cloud  and  the  radial  parabolic  density  profile  were  clear  evidence  for 
the  presence  of  a  Bose  condensed  cloud. 

Locating  the  resonances 

When  the  magnetic  field  is  swept  across  a  Feshbach  resonance  one 
would  expect  to  lose  a  condensate  due  to  an  enhanced  rate  of 
inelastic  collisions  (caused  either  by  the  collapse  in  the  region  of 
negative  scattering  length  or  by  an  enhanced  rate  coefficient  for 
inelastic  collisions).  This  allowed  us  to  implement  a  simple  pro¬ 
cedure  to  locate  the  resonances:  we  first  extended  the  field  ramp 
until  the  atoms  were  lost  and  then  used  successively  narrower  field 
intervals  to  localize  the  loss.  This  procedure  converged  much  faster 
than  a  point-by-point  search.  As  we  could  take  many  non-destructive 
phase-contrast  images  during  the  magnetic  field  ramp,  the  sharp 
onset  of  trap  loss  at  the  resonance  was  easily  monitored  (Fig.  Ic,  d)  .  | 

The  most  robust  performance  was  obtained  by  operating  the  i 
optical  dipole  trap  at  10  mW  laser  power  focused  to  a  beam  waist  of  i 
6  fxm,  resulting  in  tight  confinement  of  the  condensate  and  there-  ; 
fore  rather  short  lifetimes  owing  to  three-body  recombination j 
This  required  that  the  magnetic  field  be  ramped  up  in  two  stages:  a  : 
fast  ramp  at  a  rate  of  --lOOGms”'  to  a  value  slightly  below  that  i 
expected  for  a  Feshbach  resonance,  followed  by  a  slow  ramp  at  a  rate  ! 
between  0.05  and  0.3  G  ms‘^  to  allow  for  detailed  obser\^ation.  Near  | 
907  G,  we  observed  a  dramatic  loss  of  atoms,  as  shown  in  Figs  1  c  and 
2a.  This  field  value  was  reproducible  to  better  than  0.5  G  and  had  a  ; 
calibration  uncertainty  of  ±20  G.  ! 

To  distinguish  between  an  actual  resonance  and  a  threshold  for  | 
trap  loss,  we  also  approached  the  resonance  from  above.  Fields  | 
above  the  Feshbach  resonance  were  reached  by  ramping  at  a  fast  rate  I 
of  200  G  ms“\  thus  minimizing  the  time  spent  near  the  resonance 
and  the  accompanying  losses.  The  number  of  atoms  above  the 
resonance  was  typically  three  times  smaller  than  below.  Approach¬ 
ing  the  resonance  from  above,  a  similarly  sharp  loss  phenomenon 
was  observed  about  1 G  higher  in  field  than  from  below  (Fig.  2a), 
which  roughly  agrees  with  the  predicted  width  of  the  resonance.  A 
second  resonance  was  observed  54  ±  i  G  below  the  first  one,  with 
the  observed  onset  of  trap  loss  at  least  a  factor  of  ten  sharper  than  for 
the  first.  As  the  upper  resonance  was  only  reached  by  passing 
through  the  lower  one,  some  losses  of  atoms  were  unavoidable; 
for  example,  when  the  lower  resonance  was  crossed  at  2Gms~\ 


(0  ms)  (Time)  (70  ms) 


Figure  1  Observation  of  the  Feshbach  resonance  at  907  G  using  phase-contrast 
imaging  in  an  optical  trap.  A  rapid  sequence  (100  Hz)  of  non-destructive,  in  siiu 
phase-contrast  images  of  a  trapped  cloud  (which  appears  black)  is  shown.  As  the 
magnetic  field  was  increased,  the  cloud  suddenly  disappeared  for  atoms  in  the 
[mp  =  +1>  state  (see  images  in  c),  whereas  nothing  happened  for  a  cloud  in  the 
imp  -  -  1)  state  (images  in  d).  The  height  of  the  images  is  140  jim.  A  diagram  of 
the  optical  trap  is  shown  in  a.  It  consisted  of  one  red-detuned  laser  beam 
providing  radial  confinement,  and  two  blue-detuned  laser  beams  acting  as  end- 
caps  (shown  as  ovals).  The  minimum  of  the  magnetic  field  was  slightly  offset  from 
the  centre  of  the  optical  trap.  As  a  result,  the  condensate  (shaded  area)  was 
pushed  by  the  magnetic  field  curvature  towards  one  of  the  end-caps.  The  axiai 
profile  of  the  total  potential  is  shown  In  b. 
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about  80%  of  the  atoms  were  lost.  This,  coupled  with  the  stability 
and  finite  programming  speed  of  the  power  supplies,  limited  the 
ramp  rates  to  those  given  above. 

The  observation  of  twin  resonances  separated  by  54  ±  1 G,  with 
the  weaker  one  at  lower  field,  exactly  matches  the  theoretically 
predicted  pattern  and  thus  strongly  confirms  our  interpretation.  No 
resonance  phenomena  were  observed  in  the  \mp  =  l)  state  at  any 
field  up  to  1,000  G,  in  agreement  with  theory  which  predicted 
resonances  for  this  state  only  at  much  higher  fields. 

!  Changing  the  scattering  length 

The  trap  loss  measurements  easily  located  the  Feshbach  resonances. 
To  measure  the  variation  of  the  scattering  length  around  these 
resonances,  we  determined  the  interaction  energy  of  a  trapped 
condensate.  This  was  done  by  suddenly  switching  off  the  trap, 
allowing  the  stored  interaction  energy  to  be  converted  into  the 
kinetic  energy  of  a  freely  expanding  condensate  and  measuring  it  by 
;  time-oTflight  absorption  imaging*’'’-^  The  interaction  energy  is 
j  proportional  to  the  scattering  length  and  the  average  density  of  the 
■  condensate  («}: 

i  , 

I  Ei/N  =  —a{n)  (2) 

I  m 

’  where  N  is  the  number  of  condensed  atoms  of  mass  m.  For  a  large 
,  condensate  the  kinetic  energy  in  the  trap  is  negligible  (Thomas- 
Fermi  limit),  and  Ej  is  equal  to  the  kinetic  energy  Efcof  the  freely 
expanding  condensate  Ef^/N  =  where  vvms  is  the  root- 

mean-square  velocity  of  the  atoms.  For  a  three-dimensional  har- 


Figure  2  Observation  of  the  Feshbach  resonance  at  907  G  using  time-of-flight 
absorption  imaging.a,  Number  of  atoms  in  the  condensate  versus  magnetic  field. 
Field  values  above  the  resonance  were  reached  by  quickly  crossing  the 
resonance  from  below  and  then  slowly  approaching  from  above,  b,  The 
normalized  scattering  length  a/a  vtjN  calculated  from  the  released  energy, 
together  with  the  predicted  shape  (equation  (1),  solid  line).  The  values  of  the 
magnetic  field  in  the  upper  scan  relative  to  the  lower  one  have  an  uncertainty  of 
<0.5  G. 


monic  oscillator  potential  one  finds  (n)  ^  N(Na)~^'^  (ref.  23)  (We 
note  that,  for  a  general  power-law  potential  one  obtains 

<«>  oc  N(Na)^'\  where  1/(1  -{-  E,  l/p,)).  Thus,  the  value  of  the 
scattering  length  scales  as: 


Both  Vnns  and  Ncan  be  directly  evaluated  from  absorption  images  of 
freely  expanding  condensates.  For  a  cigar-shaped  condensate  the 
free  expansion  is  predominantly  radial,  and  so  the  contribution  of 
the  axial  dimension  to  could  be  neglected.  The  quantity 
(equation  (3)),  normalized  to  unity  outside  the  resonance,  should 
be  identical  to  ala  (equation  (1)).  This  quantity  was  measured 
around  the  resonance  at  907  G  and  is  shown  in  Fig.  2b  together  with 
the  theoretical  prediction  of  a  resonance  with  width  ^  =  1  G.  The 
data  clearly  displays  the  predicted  dispersive  shape  and  shows 
evidence  for  a  variation  in  the  scattering  length  by  more  than  a 
factor  of  ten. 

We  now  discuss  the  assumptions  for  equation  (3)  and  show  that  it 
is  approximately  valid  for  our  conditions.  (1)  We  assumed  that  the 
condensate  remains  in  equilibrium  during  the  magnetic  field  ramp. 
This  is  the  case  if  the  adiabatic  condition  d/u  <  oj,  holds  for  the 
temporal  change  of  the  scattering  length‘s,  and  a  similar  condition 
for  the  loss  of  atoms  (the  w,  are  the  trapping  frequencies).  For  the 
condensate’s  fast  radial  d>mamics  X  1.5  kHz)  this  con¬ 

dition  is  fulfilled,  whereas  for  the  slower  axial  motion 
(co,  27r  X  0.1  kHz)  it  breaks  down  close  to  or  within  the  reso¬ 
nance.  In  this  case  the  density  would  approach  the  two-dimensional 
scaling  but  the  values  for  ala  (Fig.  2b)  would  differ  by  at 

most  50%.  (2)  The  second  assumption  was  a  three-dimensional 
harmonic  trap.  If  the  axial  potential  has  linear  contributions,  the 
density  scales  instead  like  N{NaT‘'^  resulting  in  at  most  a  50% 
change  for  ala.  (3)  We  assumed  that  contributions  of  collective 
excitations  to  the  released  energ)^  were  small.  Axial  striations  were 
obser\Td  in  free  expansion  for  both  =  -hi)  and  |mf  =  -  1>  ! 
atoms  (probably  created  by  the  changing  potential  during  the  fast 
magnetic  field  ramp).  However,  the  small  scatter  of  points  outside 
the  resonance  in  Fig.  2b,  which  do  not  show  any  evidence  of 
oscillations,  suggests  that  the  contribution  of  excitations  to  the 
released  energy  is  negligible.  (4)  We  assumed  a  sudden  switch-off  of 
the  trap  and  ballistic  expansion.  The  inhomogeneous  bias  field 
during  the  first  1-2  ms  of  free  expansion  accelerated  the  axial  i 
expansion,  but  had  a  negligible  effect  on  the  expansion  of  the 
condensate  in  the  radial  direction,  which  was  evaluated  for  Fig.  2b. 

None  of  the  corrections  (l)-(4)  discussed  above  affect  our 
conclusion  that  the  scattering  length  varies  dispersively  near  a 
Feshbach  resonance.  More  accurate  experiments  should  be  done 
with  a  homogeneous  bias  field.  In  addition,  an  optical  trap  with  i 
larger  volume  and  lower  density  would  preclude  the  need  to  ramp  ! 
the  field  quickly  because  three-body  recombination  would  be 
reduced. 

The  trap  losses  observed  around  the  Feshbach  resonances  merit 
further  study  as  they  might  impose  practical  limits  on  the  possi¬ 
bilities  for  varying  the  scattering  length.  An  increase  of  the  dipolar 
relaxation  rate  near  Feshbach  resonances  has  been  predicted^  ^  but  i 
for  atoms  in  the  lowest  hyperfine  state  no  such  inelastic  binary  : 
collisions  are  possible.  Therefore,  the  observed  trap  loss  is  probably  ! 
due  to  three-body  collisions.  In  this  case  the  loss  rate  is  characterized  ! 
by  the  coefficient  K^y  defined  as  N/N  =  —  K^{tr).  So  far,  there  is  no 
theoretical  work  on  near  a  Feshbach  resonance.  An  analysis  based 

on  Fig.  2  shows  that  increased  on  both  sides  of  the  resonance, 
because  the  loss  rate  increased  while  the  density  decreased  or  stayed 
constant.  In  any  case,  the  fact  that  we  observed  Feshbach  resonances 
at  high  atomic  densities  ('^lO^^cm”^)  strongly  enhanced  this  loss 
process,  which  can  be  avoided  with  a  condensate  at  lower  density  in 
a  modified  optical  trap.  Control  of  the  bias  field  with  a  precision 
better  than  10“^  will  be  necessary  to  achieve  negative  or  extremely 
large  values  of  the  scattering  length  in  a  stable  way. 
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A  tunable  condensate 

I  We  have  observed  two  Feshbach  resonances  for  Bose-Einstein 
^  condensates  of  sodium  through  the  abrupt  loss  of  atoms,  and 
I  obtained  strong  evidence  for  a  dispersive  variation  of  the  scattering 
!  length  by  a  factor  of  more  than  ten.  ‘Tuning’  of  the  scattering  length 
'  should  become  an  important  tool  for  ‘designing’  atomic  quantum 
!  gases  with  novel  properties;  for  example,  to  create  ideal  Bose- 
I  Einstein  condensates  with  nearly  zero  scattering  length,  and  to 
I  obtain  a  detailed  pirture  of  the  collapse  of  a  condensate  with 
!  negative  scattering  length,  which  is  so  far  not  fully  understood. 

!  Tuning  the  scattering  length  can  also  be  used  to  vary  interactions 
'  between  different  species”"^  and  thus  control  the  phase  diagram  of 
'  multi-component  condensates,  possibly  switching  from  interpene- 
I  trating  superfluids  to  phase  separation^^  Feshbach  resonances  may 
!  also  be  important  in  atom  optics,  for  modifying  the  atomic  inter¬ 
im  actions  in  an  atom  laser,  or  more  generally,  for  controlling  nonlinear 
-  coefficients  in  atom  optics  with  coherent  beams  of  atoms.  □ 

i 
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